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PREFACE 


This  report  has  been  prepared  by  the  Cornell  Aeronautical 
Laboratory,  Inc.  (CAL)  for  the  Office  of  Aeronautics  Programs,  Office  of  Naval 
Research,  Department  of  the  Navy.  The  results  presented  herein  represent  an 
annual  report  of  the  tasks  associated  with  Task  Number  NR  220-041-1  of  Contract 
No.  N00014-72-C-0103. 


The  Search  Director  concept  originated  at  CAL.  It  was  initiated 
as  an  internal  research  project  as  a  result  of  various  AS! v'  studies  conducted 
at  CAL  which  were  concerned  with  ASW  surveillance  systems  and  VP  aircraft 
operations  including  detection  and  localization  of  target  submarines.  This 
study  for  ONR  is  the  first  phase  of  a  continuing  effort  to: 

1.  Formulate  the  mathematical  relationships  and  computer  soft¬ 
ware  programs  that  form  the  basis  of  the  Search  Director 
concept. 

2.  Demonstrate  the  feasibility  and  applicability  of  the 
concept  for  use  by  ASW  aircraft  crews  in  detecting  and 
localizing  submarines. 

An  abstract  of  this  study  is  contained  in  DD  Form  1173  which  is  included  in  the 
ii  and  iii  pages  of  »  is  report. 

This  study  has  been  co-sponsored  by  Aeronautics  Programs, 
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Appreciation  is  expressed  to  Cimirs.  E.  R.  Doering  and  W.  L.  Smithf  Code  461  and 
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study  has  been  performed  under  the  technical  direction  of  Cmdrs.  E.  R.  Doering 
and  W.  L.  Smith. 
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I.  SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 

A.  STUDY  SUMMARY 

1.  Purnose 


This  report  summarizes  the  significant  results  achieved  in  an 
exploratory  study  to  determine  the  feasibility  and  applicability  of  a  Search 
Director  concept  for  use  by  ASW  aircraft  crews  in  determining  optimal  search 
tactics  and  deployment  of  sonobuoys  to  detect  and  localize  target  submarines. 
Imprecise  target  information  from  acoustic  sensor  surveillance  systems  create 
problems  in  successfully  conducting  these  operations.  The  purpose,  therefore, 
of  the  Search  Director  is  to: 

a.  Assemble  this  information  in  terms  of  estimated  target 
neading,  motion,  initial  size  and  orientation  of  the 
target  uncertainty  area,  time  late,  growth  of  the 
uncertainty  area,  ASW  aircraft  search  time,  etc.  into 
coherent  and  readily  computer-processed  forms, 

b.  Analyze,  process  and  update  this  information,  and 

c.  Provide  the  Tactical  Airborne  Control  and  Coordination 
Officer  (TACCO)  with  a  dynamic  probability  density  map 
display  of  the  target’s  location. 

Probability  density  functions  are  employed  as  a  measure  of  the  uncertainty  of 
the  exact  location  of  the  target  and  its  speed  and  heading.  These  probability 
density  distributions  also  provide  measures  for  estimating  the  growth  of  the 
uncertainty  area  with  time. 

The  Search  Director  principle  makes  use  of  negative  information 
that  the  target  has  not  been  detected  up  to  the  present  time  by  deployed 
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sonobuoys.  As  areas  covered  within  the  overall  target  uncertainty  area  by 
each  of  the  sonobuoys  are  searched  and  the  target  remains  undetected,  the 
probability  that  the  trrget  is  in  the  searched  area  decreases  whereas  the 
probability  that  it  is  elsewhere  in  unsearched  areas  increases.  A  dynamic 
display  of  this  probability  density  pattern  over  the  expanding  uncertainty 
area  shows  where  the  target  might  be  at  any  particular  time  (highest 
probability  density)  and  indicates  where  the  next  sonobuoys  should  be  dropped. 
The  TACCO  can  utilize  the  Search  Director  to  examine  various  trial  sonobuoy 
patterns  and  monitoring  schedules  on  a  real  or  accelerated  time  basis  and 
select  the  most  efficient  sonobuoy  patterns  prior  to  actual  deployment  of  the 
sonobuoys. 


In  order  to  accomplish  previously  stated  objectives,  it  has 
been  necessary  in  this  study  to: 

a.  Determine  the  mathematical  formulations  required  for 
developing  computer  software  programs  to  achieve  this 
capability, 

b.  Program  these  mathematical  formulations  for  use  on 
computers,  and 

c.  Demonstrate  the  feasibility  and  applicability  of  employing 
the  concept  on  a  real  and  accelerated  tine  basis  by  dynamic 
presentation  of  results  on  displays  and  hard  copy. 

2.  Approach 

a.  Analysis 


Basic  mathematical  relationships  which  underly  the  target 
uncertainty  area  probability  density  maps  are  developed.  These  include 
analytical  procedures  for  determining  initial  target  uncertainty  area 
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probability  density  maps.  Also  included  are  subsequent  growth  factors  and 
changes  in  these  maps  with  elapsed  time  for  differing  inputs  on  target  motion 
and  effects  caused  by  the  inclusion  of  various  patterns  and  numbers  of 
sonobuovs.  \  computer  program  incorporating  these  relationships  and  procedures 
is  developed  for  use  in  determining  and  displaying  exact  solutions  for  various 
tactical  situations.  A  Monte  Carlo  computer  program  is  also  developed  utilizing 
various  algorithms  formulated  from  the  basic  mathematical  relationships.  The 
purpose  of  this  program  is  to  substantially  reduce  computing  time  and  core 
storage  requirements  while  maim  lining  output  results  that  closely  approximate 
those  for  exact  solutions. 

The  feasibility  and  applicability  of  the  Search  Director 
concept  are  demonstrated  utilizing  additional  computer  programming,  the  CAL 
Flying  Spot  Scanner  and  a  large  CRT  tc  dynamically  display  the  probability 
density  map^.  \n  analog/hyhrid  computer  system,  program  and  CRT  are  also 
utilized  to  obtain  and  display  these  maps.  Additional  subroutines  are 
formulated  as  adjuncts  to  the  Monte  Carlo  computer  program  Tor  use  in  deter¬ 
mining  sub-optimal  sonobuoy  drop -point  patterns.  Several  techniques  are 
investigated  and  compared  for  displaying  the  probability  density  maps  on  a  CRT 
in  various  shades  of  gray  nr  color.  Two  computer  programs  have  been  formulated 
to  generate  the  map  displays  as  hard  copy  on  paper. 

b.  Tactical  Cases  Analyzed 


A  comparison  is  made  of  distribution  results  from  the  Monte 
Carlo  computer  program  with  those  from  the  Exact  Solution  program  for  an 
example  problem.  Three  methods  are  analyzed  and  results  are  compared  for 
inserting  a  field  of  sonobuovs  into  the  Monte  Carlo  program.  Mathematical 
analyses  are  presented  for  development  of  the  start  and  end  grids  and  the 
Probability-of-miss  (PM)  grid. 
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I 


Three  methods  are  analyzed  for  assigning  velocity  components 
to  Monte  Carlo  (MC)  submarines.  Method  I  utilizes  a  random  number  generator. 
Method  II  utilizes  equaily  spaced  headings  of  MC  submarines  between  0  and  2  Tr  , 
whereas  in  Method  III  the  headings  are  spaced  at  equal  probability  increments 
between  0  and  2  7r  .  A  comparison  is  made  of  distribution  results  from  the 
Monte  Carlo  program  utilizing  Methods  I  and  II. 

A  biased  heading  tacticai  situation  is  analyzed  in  which  the 
target  submarine’s  velocity  is  known  and  it's  heading  is  a  random  variable 
which  is  described  by  a  probability  distribution.  A  comparison  for  this  case 
is  made  utilizing  Method  II  (equally  spaced  headings)  and  Method  III  (equal 
probability  heading  increments). 

A  tactical  situation  is  analyzed  in  which  the  target  submarine's 
velocity  is  a  random  variable.  Random  number  generator  and  expected  value 
techniques  are  utilized  in  assigning  velocity  components  to  MC  submarines. 
Computer  results  are  compai*ed  for  a  Rayleigh  distribution  in  velocity  and 
uniform  heading  utilizing  the  random  number  generator  and  expected  value 
methods  with  the  exact  solution.  Additional  results  are  compared  using  the 
expected  value  method  in  the  first  case  where  the  number  of  MC  submarines 
assigned  to  each  sampled  velocity  is  small  and  in  the  second  case  where  it  is 
large. 


A  subroutine  has  been  added  to  the  Monte  Carlo  computer  program 
for  use  as  a  possible  improved  aid  in  determining  sub-optimal  sonobuoy  drop 
points.  A  correction  factor  has  been  added  to  the  program  to  overcome  ?. 
problem  caused  by  the  use  of  a  rectangular  coordinate  system.  Sample  computer 
resuHs  are  compared  without  and  with  the  use  of  this  factor. 

In  a  sensitivity  analysis,  computer  results  are  compared  for 
changes  in  the  probability  density  maps  caused  oy  variation  of  the  input 
values  of  the  principal  parameters.  These  include  the  effects  of  varying 
the  following: 
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(1)  Start  grid  size, 

(2)  Variance  in  target  submarine  heading  and  velocity, 

(3)  Number  of  MC  submarines  per  start  grid. 

(4)  Search  aircraft  time  late  in  arriving  at  datum, 

(5)  ASW  aircraft  search  time, 

(6)  Number  of  MC  submarines  assigned  per  grid  for  a  uniform 
target  heading  and  a  Rayleigh  distribution  in  velocity,  and 

(7)  Target  velocity  on  probability  of  detection. 

Also  cor. pared  are  the  effects  of  linear  smoothing  between  grid  points,  biasing 
target  submarine  heading  and  accumulated  probability  density  of  target 
position  with  search  time. 

3.  Principal  Results 

Two  computer  programs  have  been  developed  toward  implementation 
of  the  Search  Director  Concept. 

a.  E.ract  Solution  computer  program,  and 

b.  Monte  Carlo  computer  program. 

Comparison  of  Monte  Carlo  Program  With  Exact  Solution  Results 


Probability  density  distributions  obtained  utilizing  the  Monte  Carlo  method 
closely  approximate  those  obtained  from  the  exact  solution  method.  The  Monte 
Carlo  computer  program  has  considerably  less  run  time  for  specified  tactical 
situations  than  the  Exact  Solution  computer  program  (see  Section  III,  page  III-ll). 
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A  comparison  of  sample  results  from  the  exact  solution  and  the 
Monte  Carlo  computer  program  indicates  that  the  similarity  in  probability 
density  maps  appears  to  be  excellent.  The  computing  time  on  the  IBM  370/165 
computer  for  the  exact  solution  is  92  seconds  compared  to  less  than  2.5 
seconds  for  the  Monte  Carlo  program  (see  Section  IV,  page  IV-!). 


For  a  tactical  situation  in  which  the  velocity  of  a  target  is 
known  and  its*  heading  is  uniformly  distributed,  a  one-step  random  walk 
(change  in  target  heading  during  a  specified  time  interval)  most  closely 
approximates  the  exact  solution  distribution,  whereas  considerable  deviation 
occurs  for  a  2-step  and  4-step  random  walk.  On  the  other  hand,  for  a  tactical 
situation  in  which  the  target  velocity  components  have  a  normal  distribution 
with  zero  mean  and  variance  =  1,  results  for  both  a  one  and  4-step  random  walk 
closely  approximate  those  for  the  exact  solution  (see  Section  III,  page  III-17). 

Comparison  of  Methods  for  Inserting  Sensors  in  Monte  Carlo  Program 

Of  three  methods  analyzed  for  Inserting  a  field  of  sonobuoys 
into  the  Monte  Carlo  computer  program,  Method  I  provides  the  best  results. 

In  this  method,  "looks"  by  each  sonobuoy  are  monitored  at  equal  time 
increments  so  that  one  PM  grid  can  be  used  for  all  sonobuoys  and  looks.  A 
two  dimensional  grid  is  then  calculated  as  a  function  of  x  and  y  offsets  of 
the  target  from  each  sonobuoy  (see  Section  III,  page  I 11-14) . 


The  correction  factor  which  has  been  applied  to  the  subroutine 
for  determining  sub-optimal  sonobuoy  drop  point,  determination  anpears  to 
provide  results  which  compare  favorably  with  those  which  are  expected. 

However,  a  complete  analytical  justification  for  this  factor  has  not  been 
made  (see  Section  III,  page  1 1 1-40) . 
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Comparison  of  Methods  for  Assignment  of  Submarine  Velocity  Components  in 
Monte  Carlo  Program 


Of  three  methods  analyzed  for  assigning  velocity  components  to 
MC  submarines,  Method  II  is  best.  In  this  method  MC  submarine  headings  are 
equally  spaced  between  0  and  2  Tf  and  results  correspond  closely  to  those 
obtained  from  the  exact  solution.  The  results  indicate  that  equally  spaced 
headings  allow  the  minimum  number  of  MC  submarines  to  be  utilized  in  each  grid. 
In  Method  III,  the  headings  of  the  MC  submarines  are  spaced  at  equal 
probability  increments  between  0  and  2  Tf  .  Results  for  the  biased  heading 
case  indicate  that  Method  II  provides  a  better  represents ' ion  in  the  direction 
of  the  biased  heading  while  equally  spacing  the  MC  submarines  provides  a  better 
representation  in  headings  of  low  probability  (see  Section  III,  page  III-30). 


The  expected  value  method  of  assigning  velocity  components  is 
much  better  than  the  random  number  generator  method.  Results  indicate  that 
the  expected  value  method  closely  approximates  those  shown  for  the  exact 
solution,  whereas  they  deviate  considerably  utilizing  the  random  number 
generator  method  (see  Section  III,  page  III-37). 

Monte  Carlo  Program  Start  Grid  Size,  Linear  Smoothing  and  Submarine  Assignment 
Considerations 


It  is  important  to  limit  the  size  of  the  start  grid  to  the 
minimum  required.  Doubling  the  size  of  the  grid  increases  the  computer  run 
time  of  the  Monte  Carlo  program  by  more  than  a  factor  of  four.  The  results 
indicate  that  a  2  cr  bounding  appears  satisfactory  for  use  in  analyses. 

This  encompasses  approximately  90°6  of  the  probability  density  area  describing 
the  SPA  (see  Section  IV,  page  IV-1). 
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The  use  of  a  linear  smoothing  technique  can  reduce  computing 
time  by  a  factor  of  approximately  four  for  a  specified  size  grid.  Very  little 
difference  in  probability  density  is  shown  on  a  display  in  comparing  the  case 
where  smoothing  is  employed  versus  the  case  in  which  no  smoothing  is  utilized 
(see  Section  IV,  page  IV-16). 

It  is  important  to  limit  the  number  of  MC  submarines  assigned 
to  each  grid  since  computing  time  on  the  IBM  370/ 16S  computer  increases 
linearly  as  the  number  of  MC  submarines  increases.  For  a  probability  density 
of  submarines  position  with  known  constant  velocity  and  uniform  heading,  the 
computing  time  for  40  submarines/grid  is  4.1  seconds  whereas  it  is  2.7 
seconds  for  24  submarines/grid.  A  similar  case  for  a  Rayleigh  distribution  in 
velocity  and  uniform  heading  requires  6.53  and  11.64  seconds  respectively. 

These  times  are  for  the  case  in  which  the  size  of  the  start  grid  is  determined 
by  a  2  <r  boundary.  It  appears  that  from  18  to  24  MC  submarines/grid  is 
sufficient  to  provide  a  good  probability  density  map  (see  Section  IV,  page  IV-10). 

Analog/Hybrid  Computer  System 

Results  of  the  analog/hybrid  computer  system  study  to  obtain 
probability  density  maps  indicate  that  a  general  purpose  (CP)  computer  could 
be  used  in  conjunction  with  an  Analog  Computer  (a  Hybrid  system) .  This 
system  would  use  the  G.P.  Computer  primarily  for  storage  while  the  analog 
computer  would  be  used  to  compute  the  probability  density  maps.  A  capability 
exists  for  the  analog  computer  to  pass  the  output  grid  through  a  non-linear 
amplifier.  This  technique  provides  a  method  for  investigating  the  effects  of 
intensifying  the  highest  probability  density  levels  on  a  CRT  display  (see 
Section  III,  page  III-41). 

Sensitivity  Analysis  Selected  Results 


Varying  search  aircraft  time  late  in  arriving  at  datum  changes 
probability  densities  considerably  as  time  late  increases  from  3  to  5  hours  to 
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a  maximum  of  7  hours.  The  results  for  varying  ASW  aircraft  search  time 
indicate  that  the  sonobuoys  in  effect  search  out  an  increasing  portion  of  the 
SPA  as  search  time  and  the  number  of  "looks'’  by  the  sonobuoy  field  increases 
from  9  to  13  hours  and  fron  8  to  24  "looks”  respectively  (see  Section  IV, 
page  IV-12). 


For  a  probability  density  of  submarine's  position  with  random 
velocity  and  heading,  increasing  the  variance  variation  in  heading  distribution 
narrows  the  high  probability  density  area  from  a  circle  to  an  increasingly 
narrower  ellipse.  For  a  probability  density  of  submarine's  position  with 
random  velocity  and  uniform  heading,  increasing  the  variance  variation  in 
velocity  distribution  spreads  the  hig1'  probability  areas  more  toward  datum 
than  away  from  it  (see  Section  IV,  page  IV-6). 


A  comparison  of  results  of  the  effects  of  changes  in  target 
velocity  on  the  probability  of  detection  show  che  increased  capability  of  the 
sonobuo\s  to  detect  the  target  submarine  as  its'  speed  increases  beyond  9  knots. 
This  is  indicated  by  changes  in  the  probability  density  in  the  second  quadrant 
for  the  variable  PM  grid  compared  to  that  for  the  constant  PM  grid  (see 
Section  IV,  page  IV- 19) . 

The  effect  of  biasing  target  heading  is  to  shift  the  high 
probability  density  area  of  probable  target  positions  to  a  sector  located 
outside  the  ring  of  deployed  sonobuoys  (see  Section  IV,  page  IV-24). 


CONCLUSIONS 


B. 


The  following  principal  conclusions  are  made  based  on  results 
achieved  in  attaining  the  objectives  and  accomplishing  the  five  major  tasks  of 
this  study  program: 

1.  An  assessment  of  the  principal  findings  of  this  study 
indicates  that  the  Search  Director  concept  is  technically 
feasible  and  has  general  applicability  for  use  by  ASK 
personnel  in  determining  optimum  search  patterns  for 
detecting  and  localizing  submarines. 

2.  Various  simulated  searches  utilizing  sonobuoys  deployed 

in  submarine  target  uncertainty  areas  have  been  successfully 
demonstrated  on  a  real  and  accelerated  time  basis  by  visual 
display  of  the  Search  Director  program  computer  output  on 
a  CRT,  CAL  Flying  Spot  Scanner  and  hard  copy  utilizing  an 
IBM  370-model  165  computer  printer  in  a  strike-over  process. 

3.  Based  on  limited  research  conducted  during  this  study  on 
various  types  of  display  presentations  in  shades  of  gray 
or  color  it  is  concluded  that  only  a  few  levels  of  either 
gray  or  color  may  be  necessary  in  operations  required  for 
determining  optimum  search  patterns. 

4.  The  computer  printer  strike-over  process  of  printing  a 
probability  density  map  on  paper  appears  to  be  an 
attractive  and  useful  technique  since  the  displayed 
results  can  be  readily  interpreted  and  understood  by  an 
operator.  Such  a  display  is  particularly  advantageous  in 
vehicles  where  a  CRT  is  not  available. 
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5.  Successful  development  and  computer  programming  of  both 
exact  mathematical  relationships  and  simplified  algorithms 
(monte  Carlo  program  for  computer  application)  have  been 
achieved  for  use  in  determining  (1)  initial  uncertainty 
area  probability  density  distributions  and  redistributions 
due  to  growth  of  the  area  with  time  and  differing  assumptions 
of  target  motion,  (2)  redistributions  of  probability  density 
of  target  location  within  the  uncertainty  area  as  a  result 

of  no  detections  (negative  information)  from  emplaced 
sensors  and  (3)  combinations  of  (1)  and  (2). 

6.  The  small  differences  in  probability  density  distributions 
obtained  from  the  exact  solution  computer  program  versus 
those  from  the  simplified  Monte  Carlo  program  are  considered 
to  be  negligible  and  thus  completely  justify  the  use  of 

the  flonte  Carlo  program  for  the  Search  Director  concept. 

By  using  the  Monte  Carlo  program  computer  time  can  be 
reduced  by  up  to  a  factor  of  35  over  that  required  for  the 
exact  solution  computer  program. 

7.  The  present  simplified  computer  program  of  approximately 
16,000  words  could  be  reduced  to  approximately  8000  words 
by  utilizing  half-word  techniques.  Additional  reductions 
could  be  achieved  by  time  sharing  of  some  of  the  dimensional 
statements  in  the  program.  Thus,  it  is  concluded  that 

the  smaller  program  could  be  accommodated  and  would  not 
overload  ASW  computer  systems  such  as  that  in  the  P3-C 
aircraft  or  the  system  being  contemplated  for  the  P3-C 
1975  update  program. 

8.  The  use  of  negative  information  by  the  Search  Director 
concept  (i.e.,  the  fact  that  deployed  sonobuoys  have  not 
detected  a  target  up  to  the  present  time)  has  a  substantial 
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effect  on  redistribution  of  probability  densities  to 
unscarched  portions  of  an  uncertainty  area  and,  therefore, 
on  where  additional  sonobuoys  should  be  emplaced  in  these 
areas . 

Maximum  use  should  be  made  of  linear  smoothing  techniques 
betiveen  calculated  points  in  the  Monte  Carlo  computer 
program  because  (11  computing  time  can  be  substantially 
reduced  by  up  to  a  factor  of  four  and  (2)  very  little  if 
any  differences  (with  or  without  smoothing)  can  be  seen 
on  probability  density  map  displays. 

If  computer  run  time  is  a  critical  consideration  in  actual 
ASW  operations,  it  is  important  to  limit  start  grid  size 
in  the  Monte  Carlo  program  to  only  that  required  for  the 
tactical  situation  at  hand  since  a  doubling  of  grid  size, 
for  example,  increases  run  time  by  more  than  a  factor  of 
four. 

It  is  important  to  limit  the  number  of  Monte  Carlo  submarines 
assigned  to  each  grid  in  the  program  since  computer  time 
increases  linearly  with  increases  in  the  number  of  assigned 
submarines.  Assigned  values  of  from  18  to  24  submarine 
per  grid  are  sufficient  to  provide  high  quality  probability 
density  maps. 

In  general,  results  of  sensitivity  analyses  show  that 
varying  parameters  such  as  time  late  to  datum,  submarine 
target  speed  and  associated  radiated  noise,  acoustic  sensor 
detection  capability,  variances  in  target  heading  and 
velocity  distributions,  etc.  can  have  a  marked  effect  in 
changing  probability  density  distributions  with  time, and, 
therefore,  search  tactics. 
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C.  RECOMMENDATIONS 


The  following  recommendations  are  made  as  a  result  of  the 
principal  findings  of  this  study: 

1.  In  addition  to  the  P3-C  aircraft  and  P3-C  1975  update 
program  an  investigation  should  be  made  of  the  need  for 
and  desirability  of  utilizing  the  Search  Director  computer 
program  and  associated  display  in  S3-A  aircraft,  and  as  a 
command  and  control  tool  by  the  shore  commander,  shore- 
based  Tactical  Support  Centers  and  in  surface  ships  that 
include  CV's,  LAMPS  equipped  destroyers  and  the  Sea 
Control  Ship. 

2.  Upon  establishment  of  specific  Navy  needs  for  employing 

the  Search  Director  computer  software  program  and  associated 
display  in  ASW  operations,  an  investigation  should  be 
made  to  determine  the  compatability  of  the  program  with 
Navy  ASW  computer  systems  in  terms  of  computer  core 
storage  limits,  format,  computer  language  etc. 

3.  Utilizing  results  of  the  above  investigation  as  a  basis, 
conduct  analyses  to  make  required  (1)  improvements  and 
extensions  to  the  Search  Director  computer  program,  (2) 
reductions  in  program  size  and  (3)  revisions  to  program 
in  terms  of  computer  language,  format  etc. 

4.  Formulate  a  test  program  and  conduct  an  at-sea  experimental 
evaluation  of  a  prototype  Search  Director  software  computer 
program  and  display  during  ASW  exercises. 

5.  During  the  interim,  effort  should  be  continued  to  (1)  improve 
and  expand  the  Search  Director  computer  program  and  (2) 
determine  a  preferred  method(s)  for  displaying  results  to 

an  operator. 


1-13 


II. 


INTRODUCTION' 


A.  BACKGROUND 

In  certain  U.  S.  Navy  ASW  missions,  a  problem  arises  in  search¬ 
ing  for  and  locating  a  submarine  target  utilizing  imprecise  information  concern¬ 
ing  its  position,  speed,  heading,  etc.  which  is  obtained  initially  from 
surveillance  sensor  systems.  In  order  to  increase  the  probability  of  success 
in  search,  detection  and  localization  operations,  some  means  are  required  to: 

1.  Assemble  all  available  information  (such  as  target  submarine 
heading,  motion,  position  fix  area  of  uncertainty  and  orientation, 
time  late,  growth  of  area  of  uncertainty,  search  time,  etc.) 

into  coherent  and  easily  computer-processed  forms, 

2.  Update  the  processed  information  as  search  and  localization 
operations  progress,  and 

3.  Produce  a  pictorial  representation  of  current  and  futi*-  e 
status  of  search  and  localization  operations  associated  with 
various  tactical  options  which  would  be  of  assistance  in 
conducting  ASW  missions. 

The  above  situation  occurs  frequently  in  ASW  search  and  localization,  search 
and  rescue,  and  other  missions. 

In  attempting  to  accomplish  these  objectives  a  concept  called 
the  "Search  Director"  has  evolved  as  the  result  of  a  number  of  ASW  studies 
conducted  at  the  Cornell  Aeronautical  Laboratory  in  which  search  and 
localization  operations  have  been  analyzed.  This  report  presents  results 
accomplished  to  date  from  an  ongoing  ONR  study  program  that  addresses  the 
problems  of  performing  these  operations  by  an  ASW  aircraft  in  tactical 
situations  where  imprecise  and  limited  information  is  available  on  a 
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submarine  target's  location  and  motion  within  a  specified  initial  uncertainty 
area.  It  is  expected  that  improved  usage  of  all  available  information 
provided  to  ASW  aircraft  would  increase  the  probability  of  successful  search 
operations  in  terms  of  initial  detection  of  a  submarine  by  sonobuoys  which 
are  dropped  in  a  prescribed  pattern  in  an  uncertainty  area.  The  sonobuoys 
drop  points  would  be  determined  on  the  basis  of  information  displayed  by  the 
Search  Director. 

The  operations  performed  by  an  ASW  aircraft  as  assisted  by  a 
Search  Director  are  depicted  graphically  in  Figure  2-1.  The  purpose  of  the 
Search  Director  is  to  assemble  all  the  available  information  (such  as  initial 
area  of  uncertainty  (SPA),  target  speed  and  heading,  the  aircraft's  time  late 
and  tne  search  time,  etc.)  into  a  coherent  form,  analyze  this  information, 
and  pirvide  the  TACCO  with  a  probability  density  map  display  of  the  submarine's 
location.  Probability  density  functions  are  used  as  a  measure  of  the 
uncertainty  of  the  true  location  of  the  target,  and  its  speed  and  heading. 

These  distributions  provide  estimated  growth  factors  of  the  uncertainty  areas 
with  time.  At  each  stage  of  the  search  operation,  the  associated  man  display 
is  indicated  as  shown  on  Figure  2-1.  Also  shown  are  the  target  parameter 
inputs  which  determine  the  characteristics  of  the  map  displays  associated  with 
particular  tactical  situations. 

The  probability  maps  arc  displayed  on  a  video  display  screen.  Changes 
in  the  search  situation  due  to  the  passage  of  time,  the  introduction  of  aircraft 
sensors,  and  changes  in  target  motion  cause  probability  density  levels  (light 
intensities)  to  change  on  the  displayed  maps.  The  maps  may  be  projected  into  future 
time  to  further  assist  the  TACCO  in  conducting  an  effective  search. 

In  estimating  the  probability  density  map  of  the  submarine 
target's  location,  the  Search  Director  principle  makes  use  of  negative  infor¬ 
mation  that  the  target  has  not  been  detected  up  to  the  present  time  by  any  of 
the  deployed  sonobuoys.  As  areas  covered  by  each  of  the  sonobuoys  are  searched 
and  the  target  is  undetected,  the  probability  that  the  target  is  in  the 
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Figure  2-1  OPERATIONAL  EMPLOYMENT  OF  SEARCH  DIRECTOR 
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searched  area  decreases  and  the  probability  that  the  target  is  elsewhere  in 
unsearched  areas  increases.  A  dyn  mic  display  of  this  changing  probability 
density  pattern  over  the  expanding  area  of  uncertainty  shows  the  areas  of 
highest  probability  where  the  submarine  might  be  at  any  particular  time  and, 
therefore,  indicates  where  the  next  sonobuoys  should  be  dropped.  The  TACCO  can, 
therefore,  examine  vario.o  trial  sonobuoy  patterns  and  monitoring  schedules 
on  an  accelerated  time  :  ^is  and  select  the  most  efficient  sonobuoy  patterns 
prior  to  actual  deployment  of  the  sonobuoys. 

In  summary,  it  is  believed  that  a  Search  Director  could  be 
effectively  used  to  assist  tne  TACCO  of  an  ASW  aircraft  in  determining  optimum 
search  patterns  to  be  used  in  detecting  and  localizing  target  submarines.  A 
probabilistic  process  utilizing  Bayesian  statistics  would  be  employed  to 
compute  and  display  on  a  real-time  or  accelerated  time  basis  those  areas 
which  are  most  likely  to  contain  the  target.  This  information  would  assist 
the  TACCO  in  formulating  logical  tactics  and  search  patterns  in  attempting  to 
detect  and  localize  a  submarine.  Continuous  updating  of  the  displayed  infor¬ 
mation  would  permit  (1)  real-time  evaluation  of  the  effect  of  selected  tactics 
and  search  patterns  on  detection  and  localization  operations  and  (2)  changes 
in  tactics  and  search  patterns  accordingly  as  results  are  displayed. 

Thus,  a  Search  Director  could  provide  a  continuously  updated 
visual  display  of  those  areas  having  the  highest  probability  of  containing 
the  target.  The  TACCO  can  observe  these  changing  probability  densities  as  they 
develop,  and  modify  his  search  pattern  (and/or  the  operational  mode  of  his 
sensors)  to  search  out  the  most  fruitful  areas.  Once  having  detected  a 
target  Search  Director  operations  would  cease  and  localization  operations 
would  then  be  initiated. 

Application  of  the  Search  Director  as  presented  in  this  report 
has  centered  around  its  potential  use  in  P3-C  ASW  aircraft  or  an  updated  1075 
version  of  the  P3-C.  This  selection  has  been  made  primarily  for  illustrative 
purposes  in  order  to  provide  a  means  for  detailed  discussion  of  how  it  could 
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be  employed  in  ASW  search  operations.  Other  potential  applications  include 
its  use  by  S3-A  ASW  aircraft  and  as  a  command  and  control  tool  by  Navy  ASW 
shore  commanders,  ASW  tactical  support  centers  (TSC)  or  at  sea  in  surface  ships 
such  as  ASW  support  carriers  (CV's),  LAMPS  equipped  destroyers  and  the  Sea 
Control  ship. 


The  Search  Director  could  be  used  in  S3-A's  in  a  manner  similar 
to  that  described  in  this  report  for  P3-C's.  On  the  other  hand,  information 
from  a  Search  Director  located  in  the  ASW  command  and  control  center  on  a  CV 
could  be  transmitted  to  S3-A's  and  displayed  for  use  by  crews  in  the  aircraft. 

In  a  manner  similar  to  the  latter  case  Search  Director  information  could  be 
used  by  helicopters  deployed  from  destroyers  or  Sea  Control  ships  to  carry 
out  ASW  search  operations. 

The  Search  Director  could  also  be  utilized  by  ASW  shore 
commanders,  TSC's,  and  ship  command  and  control  centers  in  planning  ASW  search 
operations,  allocation  of  available  forces,  and  command,  control  and  coordination 
of  single  and  multiple  vehicle  operations.  In  other  applications  the  Search 
Director  could  be  employed  in  TSC's  to  provide  information  that  would  be 
transmitted  to  P3-C's.  In  addition,  it  could  also  be  used  as  a  tool  at  TSC's 
in  post-operational  analysis  of  ASW  missions  and  in  training  crews  in  ASW 
search  tactics  and  operations. 
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B.  STUDY  OBJECTIVES  AND  SCOPE 


The  principal  objective  of  this  study  is  to  conduct  an  exploratory 
investigation  of  the  feasibility  and  applicability  of  the  Search  Director 
concept  for  use  by  U.S.  Navy  ASW  personnel  in  determining  optimal  search  tactics 
and  sonobuoy  drop  patterns  for  detection  and  initial  coarse  localization  of 
target  submarines.  In  order  to  accomplish  this  objective  it  is  necessary  to 
(1)  determine  the  mathematical  formulations  that  are  required  for  computer 
software  programs  to  achieve  this  capability,  (2)  program  the  mathematical 
formulations  for  use  on  computers  and  (3)  demonstrate  the  feasibility  and 
applicability  of  employing  the  concept  on  a  real  and  accelerated  time  basis 
by  dynamic  presentation  of  results  on  displays  and  hard  copy. 

The  scope  of  this  study  and  analytical  approach  are  to 
determine  the  basic  mathematical  relationships  which  underly  the  target 
uncertainty  area  probability  density  mans.  These  include  analytical 
procedures  for  determining  initial  target  uncertainty  area  probability  density 
maps,  subsequent  growth  and  changes  in  these  maps  with  elapsed  time  for 
differing  assumptions  of  target  motion,  and  effects  caused  by  inclusion  of 
various  patterns  and  numbers  of  sonobuoys.  A  computer  program  incorporating 
these  relationships  and  procedures  developed  for  use  in  determining  and 
displaying  "exact"  solutions  for  various  ASW  tactical  situations. 

In  a  subsequent  effort  the  above  relationships  and  procedures 
are  modified  and/or  simplified  for  use  in  the  development  of  a  Monte  Carlo 
computer  program.  The  purpose  in  developing  this  program  is  to  substantially 
reduce  computing  time  and  core  storage  requirements  while  maintaining  output 
results  that  closely  approximate  those  obtained  from  the  Exact  Solution 
computer  program.  Sample  results  for  identical  ASW  tactical  situations  and 
associated  variations  in  the  values  of  key  input  parameters  are  compared 
utilizing  these  programs. 
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The  feasibility  and  applicability  of  the  Search  Director 
concept  are  demonstrated  utilizing  additional  computer  programming,  the  CAL 
Flying  Spot  Scanner  and  a  large  CRT  to  dynamically  display  the  target 
uncertainty  area  probability  density  maps.  Limited  efforts  are  also  conducted 
to  (1)  determiric  the  feasibility  of  this  concept  using  an  analog/hybrid 
computer  system  and  large  CRT  to  obtain  and  display  these  maps  and  (2) 
formulate  subroutines  to  the  Monte  Carlo  computer  program  for  use  as  an  aid 
to  the  TACCO  in  determining  sub-optimal  ,onobuoy  drop  point  patterns. 

Several  techniques  are  investigated  and  compared  for  displaying 
the  probability  density  maps  on  a  CRT.  In  one  technique  the  amount  of  fill-in 
is  varied  between  calculated  points  in  which  the  probability  levels  are 
indicated  by  both  the  intensity  of  the  points  and  their  size.  In  a  second 
technique  the  number  of  shades  of  gray  presented  on  a  display  are  varied  from 
a  few  levels  up  to  a  maximum  of  61  levels  on  the  gray  scale  where  each  shade 
is  coded  to  represent  a  specified  probability  level  (interval).  A  third 
technique  examined  is  similar  to  the  second  technique  except  that  the  various 
colors  are  used  to  represent  the  probability  levels. 

Two  computer  programs  are  also  formulated  to  generate 
probability  density  map  displays  as  hard  copy  on  paper.  One  program  utilizes 
a  computer  printer  in  an  over-striking  process  of  various  letters  and  symbols 
to  print  out  the  maps.  A  second  computer  program  utilizes  a  Houston  plotter 
to  draw  the  maps  as  a  series  of  contour  lines  representing  the  vaiious  iso¬ 
probability  density  levels. 


1 1-7 


III. 


STUDY  ANALYSES  AND  RESULTS 


A.  BASIC  MATHEMATICAL  RELATIONSHIPS  UNDERLYING  TARGET  UNCERTAINTY 
AREA  PROBABILITY  DENSITY  MAPS 


1.  Initial  Target  lincertainty  Area  Probability  Density  Map 

In  airborne  search  and  localization  operations  the  problem 
arises  of  locating  a  submarine  target  within  an  area  of  uncertainty  such  as 
provided  by  the  Sound  Surveillance  System  (SOSUS).  This  area  is  known  as  a 
SOSUS  probability  area  (SPA).  One  approach  in  analyzing  this  problem  is  to 
interpret  the  uncertainties  which  arise  in  these  tactical  situations  in 
probabalistic  terms.  This  approach  suggests  the  development  of  uncertainty 
area  probability  density  maps  to  indicate  those  portions  of  the  total  area 
which  have  the  highest  probability  of  containing  the  target.  The  first  step 
in  this  development  is  to  determine  the  basic  mathematical  relationships  in 
terms  of  probability  functions  which  underly  these  maps.  In  developing  these 
relationships  consideration  must  be  given  to: 

a.  Probability  density  map  of  the  initial  uncertainty  area 

b.  Probability  density  maps  of  the  growth  of  the  initial 
uncertainty  area  with  elapsed  time  as  functions  of  time 
late  to  datum,  target  motion  and  the  deployment  of  acoustic 
sensors  within  the  area  by  the  ASW  aircraft. 

In  developing  the  probability  density  map  of  an  initial  target 
uncertainty  area  A,  it's  size,  shape,  orientation,  location  and  datum  must  be 
known.  In  addition,  the  kind  of  probability  density  function  which  best 
describes  the  probable  location  within  the  area,  e.g.,  f-(x,y )  must  be 
specified.  It  should  be  noted  that  the  function  f(x,y)  d-z  cLy  is  the 
conditional  probability  that  the  target  is  located  at  coordinates  X ,  y 
within  the  area  A,  given  that  the  target  is  located  within  the  area  A. 
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Having  this  information  the  incremental  probabilities  can  be  computed  for  use 
in  developing  an  initial  probability  density  map. 

2.  Effect  of  Target  Motion  on  Probability  Density  Map  Versus  Time 

Subject  to  the  constraints  on  target  motion  imposed  by  a 
particular  tactical  situation,  probability  densities  are  chosen  which  are 
reasonable  approximations  of  the  target's  probable  actual  motion.  In  develop¬ 
ing  these  mathematical  relationships  it  is  assumed  that  submarine  target 
velocity  and  heading  remain  constant  during  the  search  operation  unless  other¬ 
wise  specified.  Let  ot(V)'±Y,g2(&)A<j>  denote  the  probabilities  that  the 
target's  velocity  has  values  in  ( V }  V  +  a v)  and  it's  track  heading  has  values 
in  ($,  $  +  A<f>)  . 


The  first  step  is  to  compute  Pf(u.,v)AuAv-  the  probability  that 
the  target  is  ir.  +Au.)  *(r)v-+&v)aX  time  t  .  Referring  to  figure  3-1,  then 
Pt  is  given  by: 

Pt(u.,v-)A*Av  =  X  P  [target  moves  )  to  ( ic,  v~  )  in  time  2"  ] 

(t,y)  e  AtM 

*P  [target  initially  at  (  X,  y  )] 

where  the  summation  is  made  over  all  points  (  x,  y  )  from  which  transit  to 
(  u,  v  )  in  time  t  is  possible  as  denoted  by  {  u,  ir  ) . 

(*)  Pt(u,V')AuAir  =  4- $  (/>/tJ  •  f(z,y)  A*  Ay  Au  Air 

where 

Vt  =  p  =  \(u-z)2+(r-yf  \  1 ,  <f>  =  tan’1  (^) 

This  representation  is  valid  for  any  time  t  ,  subject  to  the  same  initial 
distributions  on  submarine  target  location,  track  heading  and  speed  described 
above. 
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3. 


Acoustic  Sensor  Employment  Effects  on  Probability  Density  Map 


The  next  step  is  to  compute  the  probability  density  map  for 
time  t2  (t2  >  tj)  given  that  an  acoustic  sensor  has  been  placed  at  ( x,  y  ) 
at  time  t1  as  shown  in  Figure  3-2  and  has  failed  to  detect  the  target. 

Let  Ct{R)  be  the  probability  of  detection  by  the  sensor  as  a 
function  of  range.  No  deviation  in  the  target's  motion  is  assumed  from  time 
t-0  to  time  t=  tz  .  If  the  target's  location  at  t-0  is  (x,  5  )  and  at  t~tz 
it  is  (  u,  v  )  then  the  target's  location  at  t-t,  is  given  by: 

(a.',  ix')  =  (x  *■  -j-te-x), 

This  means  that  the  target  is  just  the  fraction  t. / 1 z  of  the  way  between  (x,  j  ) 
and  {(t.-ir).  Thus,  the  probability  of  the  target  being  at  («,  ix  )  at  time  tz 
given  that  no  detection  has  taken  place  at  time  {■  is  given  by: 

Pt  [u,tr  |  no  detection  at  t1  ] 

=  P  [target  is  at  (a,  v  )  at  time  tz  with  no  detection  at  time  t1  ] 

4-  P  [no  detection  at  time  ] 

=  ^P[target  moves  from  {X,y  )  to  (  u,  v  )  in  time  tz  ] 

At(u,v) 

*p  [target  not  detected  at  time  t1  while  at  (  u,' ,  v*'  )  and  in 
transit  from  )  to  {u,  ir  )] 

*P  [target  at  (  x,  y  )  initially] 
divided  by 


rn 

/  .  P[target  not  detected  at  time  t1 

all  (Xu) 

at  (  x,  ij  )  at  time  t1  ] 


given  that  the  target  is 


[target  is  at  {  X,  <j  )  at  time  J 
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3-2  EFFECT  OF  DEPLOYING  AN  ACOUSTIC  SENSOR  ON  PROBABI LITY  DENSITY 
MAP  GEOMETRY 


It  should  be  noted  that  the  division  shown  in  the  equation  above  is  a  constant 
which  is  independent  of  u^-v.  Therefore,  it  need  not  be  calculated  since  its 
effect  is  to  normalize  the  total  area  to  unity. 

If  Pt^{t,y)  Az  Ay  is  the  probability  density  map  for  time 

r1  ,  then 

<**>  Pt  (u,  v)A«Av  =  E  T  9i  (P/tz)  9,  ($)  {/“  4-*(*,y)Az  AyAuAzr 

At£UtlO  Z  ' 

*  E  0-6(R))Pti(z,y)AxAy  , 

aJI(X,y>  ’ 

where  A  t2(u,v),p,<p  are  defined  as  indicated  previously,  and  R  denotes  the 
range  from  the  target  to  the  sensor.  In  the  numerator,  the  range  R  is 
measured  from  the  sensor  at  (z,  y  )s  to  the  point  {ex',  v'  )  described  above. 
In  the  denominator,  the  range  is  measured  from  the  sensor  at  (X,y  )s  to  the 
points  (  x,  y  )  over  which  the  summation  is  made. 

The  diagram  shown  previously  in  Figure  3-2  depicts  the 
probability  density  map  geometry  at  time  tz  when  the  sensor  is  deployed  at 
location  (x,  y  )«.  at  time  t,  . 

If  several  acoustic  sensors  located  at  ( z,  y  )  ,  .  .  ., 

(.Z,  y  )  ,  respectively;  report  failures  to  detect  the  target  at  time  t1  , 

then  ( 1- Cr(R )  )  is  replaced  by: 

n 

TT  0-  G  (R  ))  in  the  above  expressions, 

Kz)  '  * / 

where  R K  is  the  range  from  the  K'th  sensor  to  the  point  (  v-’  )  or  ( y  ) , 

as  described  above. 

a.  Sensor  Nondetection  at  Successive  Time  Increments 


If  an  acoustic  sensor  fails  to  detect  the  target  at  several 
times  past,  ti  ,  .  .  .,  tm  with  the  sensor  placed  at  {x  ,  y  )J  ,  .  .  ., 

{X,  9  ,  respectively,  the  expression  for  the  uncertainty  map  probability 

Pt  (  u.,  ir  )  A  u  A  v  is  modified  by  replacing  (f-  Q  (R  ))  with  ft  (l-Q(Rl<))in  the 

C  W  ir  -  i  v 
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rr. 

numerator,  and  ^  (  i-G(  R  ))  Ptf  (  z,  u  )  A  x  A  y  with  FI  (  (  RK  ))  x 

Pt  (x,  y  )  £i  Ay  )  in  the  denominator  of  (**).  As  stated  before,  R_ * 
in  the  numerator  is  the  distance  from  the  sensor  at  time  tK  to  the  point  tK/t 
of  the  way  between  (  %,  y  )  and  v).  In  the  denominator,  RK  is  the 
distance  from  the  sensor  to  the  point  (  X,  y  ) .  Pt denotes  the  probability 
density  map  generated  for  time  t ^  ,K-1  ,  .  .  . ,  m  . 

4.  Target  Motion  Changes  at  t 

The  probability  densities  associated  with  target  motion,  ^(V) 
and  g z  (  4>  )  are  replaced  by  g*  (  V  ),  g*  (  <$>  )  and  f  (  X,  y  )  Ax  A  y  is 
replaced  b  yP^{x,y  )£xAy  in  (  *  ) ,  with  t*  denoting  the  new  time  datum. 

5.  Target  Location  Probability  Density  Versus  Time 


Target  locution  probability  density  gt  {u,  v~)  at  tine  t  is 
derived  as  follows:  Probability  densities  g1  (V  ) ,  y2{  <p  )  are  given  on 
target  velocity  and  track  heading.  Mso  given  is  an  initial  target  location 
density  f  (  z,  y  ).  The  first  step  is  to  find  the  conditional  probability 


(p/O 


density  of  u,v  given  £,y.  Letting  p  =  Vt  ,  the  density  of  p  is  — 
measuring  p  and  cj>  with  respect  to  the  point  X  ,  y  ,  one  obtains: 

(P/t) 


By 


9t  (P’  9,  — —  9 z  &) 


By  referring  again  to  Figure  3-2,  it  is  evident  that  the  follow¬ 
ing  relationships  exist: 

u  =  x  +  p  cos 
v  -  y  -f-  p  s  m  $z5 


Thus , 


where 


Here, 


St  («,r,z,y)=  j-  gf(Vt) 


d(/>,  d>) 

d  (u.,  v) 


<fi  =  tan 


B(p,4>)  _  A  /d(u,  v)  _ 
3(u,ir )  /  d(p,  <f>)  '  ,P 


— - - —  is  the  .Jacobian 
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of  />  ,  4>  with  respect  to  u,v  (z,y  are  considered  to  have  fixed  values  in 
computing  the  conditional  probability  density). 


Hence, 


st  (u.vV. s)  =  FSi(f  \Lu -*)*+<*-&*} *) 3 2 

st  fa,  V)  =JJ gt  (u,  jr  |  *,  y)  f  (z,  y)  dxety 


Finally, 


(fe-v)2  Hv-tf)  z~j 


V* 


9tr  fa^  1  z>y)  t  fa  <j)d-xdy  , 

J4*-v> 

where  Ai  {u,v-  3  is  the  set  of  points  {t,y  )  such  that  transit  to  ( u,ir )  is 
possible  under  the  constraints  of  g  and  gz  . 


The  equation  for  p*  («,*■  )  Au  Air  is  an  approximation  to 
g  (a,  v)  du.  ctv  . 


B.  DESCRIPTION  OF  "EXACT"  SOLUTION  COMPUTER  PROGRAM 


1.  Method  of  Determining  Probability  Density  Man  of  Target 
Position  With  Time. 

A  computer  program  has  been  developed  which  closely  models  the 
exact  mathematical  relationships  governing  the  Search  Director  concept.  The 
principal  anproximation  in  this  program  is  that  the  continuous  probability 
densities  are  approximated  by  discrete  functions.  The  program  first  divides 
the  initial  probability  density  area  or  SPA  into  NxN  grids  and  then  calculates 
the  probability  of  the  target  being  in  each  Ax  by  Ay  grid  at  time  equal  to 
zero  when  the  target  is  initially  detected  by  SOS US  and  a  SPA  is  formed.  This 
probability  is  given  by: 

Pj  --  Ay 

where  PLJ  is  the  probability  of  the  target  submarine  starting  in  the  i  ,  j 
grid  and  -f  (  X ,  y)  is  described  by  a  joint  circular  normal  probability 
density  (SPA).  The  point  XL  ,  yt-  is  at  the  center  of  a  grid  of  size  Ax 
by  A  y  . 


m-s 


From  the  estimates  of  the  heading  and  velocity  probability 
densities,  four  grids  (V  g  ,  B- 1,2,3  or  ‘j-  )  are  calculated  for  the  probability  of 
the  target  starting  in  grid  L  , j  and  going  to  grid  m,  n  in  time  t  .  The  four 
grids  are  calculated  for  the  following  constraints  on  Lt  j  and m,  r?  : 


Grid  1 

Grid  2 

Grid  3 

Grid  4 


m  i  i 

n  -  j 

m  -  L 
n  ±  j 

m  ■£■  i 

n  >  j 

M  -c.  i 

n  *  j 


2.  Method  of  Inserting  Sonobuoy  Field 


The  sonobuoy  field  is  inserted  into  the  computer  model  by  first 
calculating  a  composite  probability  of  miss  grid  for  all  the  sonobuoys  at  each 
"look"  interval.  The  position  ( -ktJL)  of  the  target  going  from  grid  ifj  to  m } 
n  at  the  time  of  the  ath  look  is  given  by: 

=  i  +  ~  (m  -  C) 

*  a  j  *  y-(n'j) 

Where  is  the  time  of  the  ath  look  and  T  is  the  time  for  which  the 
probability  density  is  to  be  calculated.  The  above  equations  are  valid  for 
the  assumption  that  a  covert  search  operation  is  conducted  and  thus,  target 
heading  and  speed  remain  constant.  The  composite  probability  of  miss  for  a 
target  starting  in  grid  i)  j  and  going  to  gridm;n  is  given  by: 

A ILOOK 

pm  -TT  zj 

0.-1 

where  NLOOK  is  equal  to  the  number  of  discrete  "looks"  by  the  sonobuoy  field 
and  is  the  composite  probability  of  miss  grid  for  the  ath  look. 
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To  find  the  probability  of  the  target  submarine  being  in  the 
m,n  grid  at  the  update  time  T  ,  the  following  operation  is  performed. 


where Vg  is  the  grid  representing  the  velocity  and  heading  distributions. 
C.  DESCRIPTION  OF  MONTE  CARLO  COMPUTER  PROGRAM 


1.  Method  of  Determining  Probability  Density  Map  of  Target 
Position  With  Time 


Algorithms  have  been  developed  to  "grow"  the  initial  probability 
functions  such  that  a  general  purpose  computer  can  be  used  to  provide 
probability  density  maps  of  target  uncertainty  areas  on  an  accelerated  time 
basis.  One  of  these  algorithms  has  been  developed  to  grow  the  initial 
probability  grid  by  a  Monte  Carlo  (MC)  method.  The  computer  program  solves 
the  initial  probability  density  distribution  (SPA)  in  the  same  manner  described 
previously  for  the  exact  solution.  The  initial  distribution  is  assumed  to  be 
a  joint  normal  distribution  in  which  the  X  and  y  components  are  assumed  to 
be  independent  but  their  variances  need  not  be  equal. 

A  number  of  Monte  Carlo  elements  (MC  submarines)  are  then  placed 
in  each  Ax  by  Ay  grid.  The  same  number  of  MC  submarines  is  inserted  in  each 
Ax  by  Ay  grid.  From  the  velocity  and  heading  distribution,  x  and  y  components 
of  velocity  are  assigned  to  the  MC  submarines.  The  same  set  of  velocity 
components  are  used  for  the  MC  submarines  in  each  Ax  by  A  y  grid.  A 
discussion  of  the  method  used  in  this  model  to  obtain  the  velocity  components 
is  presented  in  Subsection  C-6  of  Section  III  of  this  report. 

The  computer  program  then  converts  the  velocity  components  into 
the  number  of  Ax  by  Ay  grids  in  which  the  MC  submarines  move  at  the  update 
time  T  .  Knowing  the  MC  subs  starting  grid  and  the  number  of  grids  that  the  MC 
submarines  travel  in  the  x  and  y  direction,  the  coordinates  of  the  MC  submarines 
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end  grid  at  the  update  time  T  may  he  found.  A  count  is  then  performed  on  the 
number  of  MC  submarines  in  each  grid  at  the  update  time  weighted  by  the 
probability  of  the  MC  submarine  starting  grid.  This  count  is  normalized  In- 
dividing  by  the  number  of  'HI  submarines  placed  in  each  Ax  by  Ay  start  grid. 

2.  Comparison  of  Monte  Carlo  Computer  Program  Sample  Results  With 
Hxact  Solution  For  Initial  Probability  Density  Distribution 


Figure  5-3  presents  results  from  a  computer  run  from  the  Monte 
Carlo  program  together  with  those  for  the  exa-.t  solution.  The  initial 
distribution  (SP\)  is  a  joint  circular  normal  with  a  standard  deviation  of  one 
(  <j"5Pa  =  1  )  for  the  %  and  y  components.  The  velocity  components  (  vg  ,  Vy  ) 
were  each  generated  from  a  A '{0,  1  )*  random  number  generator.  The  curve 
shown  for  the  Monte  Carlo  method  has  been  calculated  for  the  time  t-1  and  is 
a  cross  section  of  the  joint  distribution  passing  approximately  through  the 
mean  of  the  distribution. 


The  curve  shown  for  the  exact  solution  is  calculated  in  the 
following  manner: 


Let  cc  -  i  +  irxt 
v=  y  +  Vyt 


ilence,  ( a,  v )  gives  the  position  of  the  target  starting  at  y  with  velocity 
components  ( )  at  the  time  t  . 


The  random  variables  x  and  y  are  described  by  the  distribution  governing  the 
S!'-\  N  (0,1  ;  0}1  )*  and  irx>  are  also  described  by  a  joint  normal. 

The  expected  value  of  u  (e  J  is  giver,  by 
£{u}  *  E  [tv-  »r  t} 

=  0 

*The  notation  N  (0>l  )  denotes  a  normal  distribution  having  0  mean  .<nd  variance 
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DISTRIBUTION  OF  A  MOVING  TARGET:  X  AND  Y  VELOCITI ES  ARE  N  (0, 1) 


The  variance  of  u.(Va.r  (u.))  is  given  by 
Va.r  |  j-  -  E  z  *■  Vx  t  | 

=  E  +  2  z  vx  t  /■  v*  t*  j- 

*  «{**}  "**{»>*} 

-  cr  -r  t  <r 

spa  tr^ 

Since  u.  is  described  by  the  sum  of  two  normal  random  variables,  u,  is  also  a 
normal  random  variable  with  a  mean  of  zero  and  a  variance  equal  to  (Ts ^  y- 1 z  . 
The  expected  value  of  iris  also  zero  and  has  a  variance  of  cSPA  + 1 2  cr ■*  . 

A  comparison  of  the  two  curves  in  Figure  3-3  indicates  that 
the  probability  density  distribution  obtained  utilizing  the  Monte  Carlo 
Method  closely  approximates  that  obtained  from  the  exact  solution.  The 
relatively  small  differences  shown  are  negligible  for  CRT  display  of  the 
probability  densities  as  variations  in  light  intensities.  Thus,  it  is 
concluded  that  the  Monte  Carlo  Method  which  has  considerably  less  run  time 
requirements  on  the  computer  is  the  preferred  technique  to  use. 

3.  Methods  For  Inserting  Sonobuoy  Field  Into  Monte  Carlo  Program 

Three  methods  have  been  investigated  to  insert  the  sonobuoy 
field  into  the  Monte  Carlo  program.  It  is  assumed  in  each  of  these  methods 
that  the  sonobuoy  field  searches  the  area  in  discrete  "looks".  This  in 
effect  assumes  that  the  output  stage  of  the  signal  processor  is  an  integrator 
and  that  this  output  is  sampled  at  the  Nyquist  rate.  Also,  it  is  assumed  that 
all  integrated  outputs  of  the  signal  processors  are  sampled  at  the  same  rate 
and  at  the  same  time. 

In  all  the  results  presented  in  this  report,  the  probability 
that  the  sonouboy  detects  the  target  at  each  of  the  discrete  "looks"  is 
assumed  to  be  of  the  form: 

P D  -  exp 
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where  R  is  the  distance  of  the  sonobuoy  to  the  position  of  the  target  sub¬ 
marine  at  the  time  of  the  ’’look"  and  oc  is  an  input  to  the  model.  This 
probability  of  detection  curve  is  not  a  limiting  factor  in  the  model  but  is 
only  used  for  convenience  in  analyses  conducted  during  this  study. 

The  concept  of  Search  Director  make  use  of  the  negative 
information  that  the  target  has  not  been  detected  which  is  given  by  1-PD 
and  is  denoted  by  Pm  . 

a.  Method  I 


In  this  method  of  inserting  the  sonobuov  field  into  the  Monte 
Carlo  program,  the  location  of  all  the  sonobuovs  and  the  "looks"  on  each 
sonobuov  that  have  been  monitored  are  put  into  machine  storage.  It  is 
assumed  that  the  "looks"  by  the  individual  sonobuovs  are  at  equal  increments 
of  time  so  that  one  PM  grid  can  be  used  for  all  sonobuoys  and  "looks".  A 
two  dimensional  grid  PM  (  i,j  )  is  then  calculated  as  a  function  of  the  z 
and  y  offsets  of  the  target  from  the  sonobuov.  A  more  detailed  description 
of  the  probability  of  miss  grid  is  presented  in  Subsection  C-S  of  Section  III 
of  this  report.  As  an  MC  submarine  moves  from  grid  C,j  to  grid  mt  n  the 
location  of  the  MC  submarine  at  each  look  by  the  sonobuoy  field  is  calculated 
from  the  following  relationship: 

ka.  =  L  * -j-  (m-i) 

By  knowing  the  position  of  the  MC  submarine  at  each  "look",  the  distance  of 
the  MC  submarine  from  each  sonobuoy  is  calculated. 

Utilizing  the  relationships  above,  the  accumulated  probability 
that  an  MC  submarine  remains  undetected  is  calculated  by  the  following 
relationship: 

fJSONQ  A/LOOK 

PM  =  TT  TT  ■pm ,  ■ 

1=1  J 
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where  pm  ■  is  the  value  taken  from  the  PM  grid  to  represent  the 
probability  that  the  ith  sonobucy  does  not  detect  the  MC  submarine  on  the 
jth  look.  Also,  N'SONO  is  the  number  of  sonobuoys  deployed  and  NLOOK  is  equal 
to  the  number  of  "looks"  by  the  sonobuoy  field.  If  the  ith  sonobuoy  is  not 
being  monitored  on  the  jth  "look",  then  pm  ..  -  1  .  Having  the  value  of 

PM  ,  the  count  of  each  MC  submarine  in  grid  rr,,  n  is  weighted  by  P  •  *  PM  . 

b.  Method  II 


Assume  that  the  initial  probability  density  distribution  of  a 
SPA  is  given  by: 


f  K  !/)  = 


2  ir  crxCy 


exp 


1  *  w 


LetVYvj0)be  the  joint  distribution  describing  the  motion  of  the  target.  By 
making  the  following  transformation: 

Wx  -  v  cos  (0) 

Vy  -  v-  s/n(6 ) 

the  joint  distribution, y  ,  Vy  ),  of  the  velocity  components  is  obtained. 

Let  the  start  grid  of  the  submarine  target  be  denoted  by  }  yj 
and  the  end  grid  by  Xm,  z n  then  the  probability  that  the  target  moves  a 
distance  %m-  in  the  z  direction  and  yn  -  y  L  in  the  y  direction  in  time  t 

is  given  by: 

S{v'x>v'y)Av-x  Av-y  -  ±  3  AzmAyn 

The  probability  of  the  submarine  starting  in  the  grid  ,  yj  is  given  by: 

&yj 

Hence,  the  probability  of  the  submarine  being  in  the  grid  2-m  ,  y n  at  time  t 
is  obtained  by  summing  over  all  possible  starting  grids: 

h(xm,yn,t)AxmAyn  =  jrZ  Z  f(Zt>yj)g(^j^,^^AZiAyjAzmAyn 

J 
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Assume  that  one  "look"  at  the  area  is  made  by  a  field  of  sonobuoys  and  that 
the  probability  that  the  sonobuoy  does  not  detect  the  target  is  given  by: 

PM  -  xa  ,  %  -y0  ) 

where  ( )  denotes  the  location  of  the  target  and  (  z0  ,  y0  )  the 
location  of  the  sonobuoy. 


The  joint  probability  of  the  target  starting  in  grid  x^,  y : 
and  going  to  grid  yn  )  and  that  the  target  is  not  detected  is  given  by: 


h  (fm  -  Vn  ,t)  At  abyn  =  -jt  Z  E  r  -y  •  Pm(xc  *  |*-  (tm  ~XL  )-Xa,  yj  +  ^L  (yn~yy)  -  ya) 


whe-'c-  rf  is  the  time  of  the  "look"  t is  the  time  of  the  update  and  K  is  a 
normalizing  constant. 


At  the  time  of  the  first  "look",  t  =  t  ,  PM  is  not  a  function 
of  the  increments  being  summed  over: 

h  (Xm>&n>  t»tf)* 

£  J 

At  the  second  look  consider: 

b  (xm,  un  ,  t--  tf  +  At)  =  PM2  Z  Z  n  (*L ,  t/i  ,t) 'g 

•  j 

The  result  of  this  method  approximates  the  motion  of  the  target  as  a  random 
walk.  That  is,  at  each  new  time  update  the  target  is  allowed  to  turn  to  a 
new  heading  governed  by  the  heading  probability  density. 

When  no  sonobuoys  are  deployed,  the  effect  of  utilizing  the 
above  method  is  that  at  the  end  of  N  steps  of  At  length,  the  variance  of  the 
resultant  probability  distribution  is  proportional  to: 

t  AtY 

i-f 
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while  the  assumption  that  the  target  maintains  a  constant  heading  and  velocity 
the  desired  variances  is  proportional  to: 

(  AtJ 

One  method  of  correcting  this  is  to  multiply  the  components  of  velocity  by  a 
constant  at  each  step  so  that  the  resultant  variance  is  equal  to  the  desired 
variance.  Thus,  a  series  of  constants,  ,  are  chosen  so  that 

//  /'J  \  2. 

Z  ia^iUT  At) 

t:1  \<-'7  / 


for  all  values  of  L 


c .  Sample  Results  Utilizing  Method  II 

Sample  results  utilizing  Method  II  are  presented  in  Figures  2-4 
and  o-5  for  the  case  in  which  sonobuoys  have  not  been  deployed.  Figure  3-4 
shows  results  for  a  tactical  situation  in  which  the  velocity  of  the  target  is 
known  and  its'  heading  is  uniformly  distributed  between  0  and  Iff.  Curves  I, 

II  and  III  show  a  comparison  of  the  probability  density  distribution  of  a 
moving  target  for  a  1-step,  2-step  and  4-step  random  walk  (changes  in  target 
heading  during  a  specified  time  interval)  respectively.  Curve  II  most  closely 
approximates  the  exact  distribution  whereas,  considerable  deviation  occurs 
for  Curves  I  and  III.  Thus,  it  is  concluded  that  a  random  walk  of  more  than 
one  step  will  alter  the  actual  probability  density  map  considerably. 

Computer  run  results  in  which  the  target  velocity  components  are 
N  (  0,  / )  are  presented  in  Figure  3-5.  Curves  I  and  II  show  a  comparison  of 
the  probability  density  of  a  moving  target  for  a  4-step  and  1-step  random  walk 
respectively  with  that  computed  for  the  exact  solution.  In  both  runs  the 
velocity  components  of  the  MC  submarines  are  multiplied  by  a  constant  at  each 
step  so  that  the  variance  is  equal  to  the  desired  variance.  On  this  basis,  it 
is  concluded  that  both  Curves  I  and  II  closely  approximate  the  curve  for  the 
exact  solution. 
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STEPS  (CURVE  I) 


AND  UNIFORM  HEADING:  RANDOM  WALK 


Figure  3-5  DISTRIBU1  !0N  OF  A  MOVING  TARGET  WHEN  X  &  Y  VELOCIT  Y  COMPONENTS 
AREN  (0,1):  RANDOM  WALK 


d.  Method  III 


This  method  approximates  the  update  probability  density  as 


follows: 


wnuno.  tt — •  r — v,  _ 

iLY. J-g-PM=  p(z,y,t)'  PM  (x,y,t) 
where  p  (z,y,  t  )  =  2  Z  f  -g  ■ 


The  method  of  obtaining  p  {X,  y ,  £  )  is  by  the  Monte  Carlo  method  described 
previously.  The  PM  grid  is  then  "grown"  separately  from  the  p  {x,y,t) 
term.  The  only  case  modeled  in  this  study  for  a  heading  distribution  of  the 
target  which  is  assumed  to  be  uniform  between  0-2  7f  . 


The  PM  grid  is  set  equal  to  the  probability  that  the  target 
is  not  detected  on  the  first  "look"  by  the  sonobuoy  field.  Between  successive 
"looks"  the  PM  grid  is  then  "grown"  by  the  following  operation: 

PM  *)  =  i  +  PM (Z,_  1 ,  yj,  t-At) 

+  yj ,  i+At)  *  PM  (X;,  yj+f, 

For  heading  distributions  other  than  for  a  uniform  heading  between  0  and  2  if  , 
the  above  operation  would  be  performed  by  weighting  each  term  in  the  sum. 


At  the  time  of  each  "look",  the  PM  grid  is  multiplied  by  the 
probability  that  the  sonobuoy  field  does  not  detect  the  target  during  that 
"look".  The  term  At  in  the  above  equation  is  set  equal  to  At/E  <|V}-  where 
is  the  grid  size  (assuming  A  X  =  Ay)  and  E  {^}is  equal  to  the  expected 
value  of  the  target’s  velocity..  Let  T  be  the  time  increment  between  "looks", 
then  the  above  operation  is  performed  N  times  between  "looks"  where  N  is 
defined  as: 

N  =  T/At 2 
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Sample  Computer  Run  Results  For  Methods  I,  II  and  III 


Figures  3-6  through  3-11  present  results  of  some  sample  computer 
runs  for  the  three  methods  of  inserting  the  sonobuoys  into  the  Monte  Carlo 
program.  The  plots  were  obtained  from  an  on-line  printer  used  in  conjunction 
with  the  IBM  370  Model  165  computer  system.  Since  the  printer  prints  10 
characters  per  inch  in  the  horizontal  direction  and  8  characters  per  inch  in 
the  vertical  direction,  some  distortion  is  present  in  these  figures.  The 
symbols  on  these  figures  indicate  various  probability  density  levels  with  the 
darker  areas  indicating  higher  probability  density  levels.  These  levels  are 
stepped  in  approximately  5%  increments. 

Each  typewriter  character  shown  on  the  figures  represents  a 
grid  of  .1  unit  by  .1  unit.  The  standard  deviation  of  the  initial  probability 
density  of  the  SPA  is  a  joint  circular  normal  with  a  standard  deviation  of  1  unit  in 
both  the  X  and  y  directions.  The  SPA  is  centered  approximately  in  the  center 
of  each  figure.  This  point  is  also  called  datum. 

The  number  of  emplaced  sonobuoys  shown  in  each  figure  is  two 
with  each  sonobuoy  having  S  ’’looks".  The  sonobuoys  are  located  1  unit  to  the 
right  and  left  of  datum.  In  all  cases,  the  target  submarine  heading  is 
uniformly  distributed  between  0  and  2  7f  . 

For  Figures  3-6,  3-7  and  3-8,  the  target  submarine  has  a  known 
velocity  of  1  unit/unit  of  time.  The  display  (elapsed)  time  is  2  units.  In 
Figure  3-6  the  sonobuoys  are  inserted  by  Method  I.  Figure  3-7  by  Method  II  and 
in  Figure  3-8  by  Method  3.  Figures  3-9  and  3-10  are  the  same  as  Figures  3-6 
through  3-8  except  that  the  velocity  component  of  the  MC  Submarines  have  a 
normal  distribution  (a/  (  0,1  )).  The  magnitude  of  the  velocity  has  a  Rayleigh 
distribution.  I.i  Figure  3-9  the  sonobuoys  are  inserted  by  Method  I,  Figure  3-10 
by  Method  II  and  in  Figure  3-11  by  Method  III,  A  comparison  of  these  and 
additional  results  indicate  that  Method  I  of  inserting  the  sonobuoys  into  the 
Monte  Carlo  program  provides  the  best  results.  For  this  reason,  a  major 
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SUBMARINE  HEADING:  UNIFORM  DISTRIBUTION  0-27 T 
SUBMARINE  VELOCITY:  1  UNIT/UNIT  OF  TIME 
ELAPSED  TIME:  2  UNITS 

NUMBER  OF  SONOBUOYS:  2  NUMBER  OF  "LOOKS":  5 
TIME  OF  "LOOKS":  1 .0, 1 .25, 1 .5, 1 .75,  2.0 


Figure  3-6  PROBABILITY  DENSITY  OF  SUBMARINE'S  POSITION  WITH  KNOWN 
CONSTANT  VELOCITY  AND  UNIFORM  HEADING:  METHOD  I 


1 1 1-22 


—  ••  ••  •  •“  -  -  •  ■  *#•»«•••• 

- - (Him  Mtm  - - 

- - =  - - 

- - ;:.r:H44M))))|l)])))||)|)tMHfti::r: - - 

- - 1  ))  }  |  )  1  1  1  |  )  1  j  i  1)  )  )  )  )  )  144444  =  :=  -= - - 

- - ::7:(li)||)||il!|lHi[||!ll!!l|||)||]UH:=:= - - 

- - -■---****  )U  )\\\\  \  h  Ui/U/ /.  Z/'/Zlllll))!)n**  =  =  =  = - - 

- - .r  =  .=+44+)Hillli//'/XXX>.X<'X/^X///,;Zllll)J)‘*  +  «- =  . - 

- .  .  .  .  =  =  _  =  4  +  4-  )  )  )j  ]]  /  ZZ  X  XX  X  A... •.A/iXXX  A  XiT  U.L  1  1  1  )  »  >  4  <  4  =  =r  =  . . . . - 

- - =  =  =  4T4  j  j  )  !  i  i/,./XAX*.r;‘4VVW«v.  ''.vy*4A/.<>  X/ZZ1  1  1  1  1  )  4  4  4  =  1=  .  .  .  . - 

- . - "  =  44*1  1  1  1  UZ  /  X  ‘  XAA  »"•:'•••  ^-rrs— or-.t-w^-.v/.AAXXXZZ/  !  1  1  1  l!  *4s".... — 

_ _ —  —  444)1)11  u /*./.,  i.iA  •  •  it=~n-;evVv.:1AAXXZ//  11)1)44=  =  =  .... — 

_ _ =  =  =  *4)))i  lz//y4.'.;.-"-trifi;=!;=f'H*  =  e3=  =  #£t£'"'AAx  xx  //  l 1)1)44  =  =  =  .  — - 

-....  =  =  =  *4 ) ) )  i  *.  .-:-==;-=s.4  *  *^7#rt:4cf+44;44fe=£r->'yy<\Axx//i  n  1)44  =  =  =  ... 

-....=-444)  )  1  1  Ll  'A  "-.-rf=S4-*.  47* £*.£****  4  tt-t  T^=5i5*4H*--',AA  X/ U  1 1  )  )  4  4  =  =  = - 

....  =  =  =  44))  11=/  .<  X  -V  4 '  • '  r.  r.  g  ?  4  *  ,*A:  i  •»  :f  £  f  £  i  £  £  £  £  f. «  £  £  t  g  F.  h  r  X  A  A  X  X  /  / 1  I  1  144  =  =  =.. 

.  .  .  .  =  =  4  4  4  )  )  i  iX/XAA  1  1  1  )  4  4  4  =  =  .. 

...-  =  =  44)111  /  t>  xA-\f===t*  >f£l£«-*5S**4i45ia*efttEEcc  '•-•AxXZZ  1 1  )  )  44=  =  =  . 
...  =  =  4  4  4  )  1  1  /;///  AA-'<---3>i'4SS41:4*«*SS«**S*4i:i1i*4>tiS;5er-"  M-Wll  1  )  )  4  4  =  =  =. 
..  =  =  =  44)  )  1 1  /  /  <  A  A  *•*  -r=;i'.t»iat*£iiSS§*8S6Eettit<£e£f.£f§Sr--’'AAXZ<’l  ,  )  )44=  =  = 
,,==itf)  J  1  /./X/.AAr--.r:4=i=?  #f4££S<S£»ifSeeEEBEf  8*41iJ%?a«=!Hn  -AAXX:/  1  )  )  4  4  =  = 

.=  =  =  4  4  4  )  1  1  /  /AXA■•yr=v==-<:4A4■ii^iieee■*eB^EEaEe^^!i^££^?•.■»1^-r•,  v  A  A  X  Z  /  1  1  )  4  4  4  =  - 
.  =  =  =  +  4)  )  ii/  xx\A  >-£-=,*  S4M.ili£SSaa8a8aBS*IiBff**ii£tffH-:4=.r->  VAAX/./  ]  !  )  )4  +  .-=  = 
.  =  =  =  4  4  )  )  1  /  /X  =  A--"--.-.->.4,f<££«*®5aBiil«iiaaaf  SS£4it«:f4)i=-jr,/,,AX  x//  1  )  >  44  =  =  = 
.  =  =444  j  1 1//  *  a  vv-‘-==fc-fft*£fisxa39saiMiiaa(J8«si4ii*>jM  =-_••»;/,*  a  x*z  11)444  =  = 
=  =  =4+ )  j  1  i/xxAA-,  'c-,=3=+H,.t£-e*s6saaBBaaaaaaaaa*£-tHi*^s*  =*  .=)  '.v.xx/  111)44-- 

=  =  =44)  )  U /XxiA:-.-,£-E  =  #»;i£«£-tS88BaB*iBSBe8BE£S*3-i3-iif===5  fcXXAX  </  11)  )44=  = 

=  =  =  44))  1 1 Z  X  X  A •' X  ,-r.=?tti^i-li*SBS*rBBaaBBfaa8aSS4S££f>.f=S=)  IhKMAAX/:  1  1)  )  44  =  = 
=  =  =  4  4  )  )  1  1/  x x A  ■:y'.<-1-,^f4:f*£ili**BBaBSBBB155BaaBJ££e£f  £4>s=i....>;yVAAX/ '  1  1  )  4  4  4  = 
=  =444 » )  i  z/xAv-i’- M-e4£i8aeaaBBBBBBBBest*i4£ft?\-k.3i=?;=‘;-(:v>*AAx.;  zi  1 1**4= 
=  =  444  )  )  1  ZZX  AA;,..;y3ji=**»  40  i«e8BlBE£BBB8fi* AX/  Z  1 1 )  )  4  4  = 
=  =  44  )  )  )  1  /xx A  A  »•  i-iH5?=4  •.•'.••'-rS**fiBa8*BBB88B'i?{=-'-X-,=J9J,4--t4snr.’'vXAAXZZl  l  )  )  4  4  = 

=  =  44  j  j  1  iz  xx aav-  1  i/-.*aaaaa3aaasa**‘-yi  1  /Affctsst—f'y  aaxzzi  1  )  >  4  *■-- 

=  =44  )  )  1  1  ZaXAA  \A  '■-,d=x=  •./  4  =  =  )Xr£CB88BBBBSB-1-=X  )  =  =  4/.-  •'.^■6  3  =?-!-!■->  AAXZZ  11  )  )  44  = 

=  =  44))  1  l/xx..  -.a*'.  1.  — =/-:iasaaattaaaa£- 1- —  =  1  sfrFs£-xvx vax  /  z  1 1 1 )  44= 

=  =44)  j  1 1  / x x /i a \ a •* x- = = -  v. .--  =  /-£aaaaaaaaaai-z= — =  1  A=)3B=-ic=xyxAx.:/ 111)44= 

=  =  44  )  )  1  1  /.  X  A  A  AM-V  ■•’-)  w  S“  '4  =  =  )X7=-£SiaBBBBBBEi=  X  )=  -  4ZyPH:!==r—  *  v  AX  /  Z  1  1 )  )  4  4  = 
=  =  44)  )  l  1ZXXAAA  -  iAZ  j  lXfctlBBBBaBBiiBlScX)  1  /  A  r  £’*■'?■  &e  *;  v  y -i  A  X  Z  Z 1  1  )  )  44  = 

=  =  44  )  )  1 1ZXX*-A  f  '-jr=.=  <;F-  -'.‘'=.«BaBBBaaBBBaE£  £  £  P  't  =  £££fc*r-!e"v'/AyZZ  1  11)44  = 
=  =  44 1 )  i  lzxx'.A  '  h*  f  4  £.tiaaaBaaaaBiBB®i.i  »-4»»-.t£:4!is-r>r:f  >•  xsx/Z  111)44  = 

=  =  4  4)  )  I  1  /  x  X  *  -If'  fit  if  £-i*  SBf  BSaiHBBBBBESEliSif  )f!=1--;- AX  /  Z  1  1  )  )  *4  = 

=  =  ♦  +  )))  1//.X-'  A’’’-  rs— -=i=s-X;  tit  l-iilS8Ba8aBiB8BBBB8B88S4£44£ A  < '  Z  1  1  )■*♦*- 
=  =  444  )  )  1ZZ  XAA  -i'<=Hgt,S>-4.4*4**!£**e«SBBaiaBBBBBa8tBJ<ifta4j4i?F-=r:c’-AA\ZZ  1  )  )  4  4  *- 
=  =  =44  )  )  1ZZXAA-1,  -,r.=.p44?i4iiS*ttaaafa«ai»*SBBEBt3-*J;ia4-4!f!'=3-r'-'.AA>;.Z  1  )  i  4  +  =  = 

—  =  =  44 1  j  1 1  / x x a *i *  rHFga.vUiiiaBBaaaiaaaaaBfaaBai^iieae^TXhrt-if-  aaxzzi  ) )  4  4  - 
=  =  =♦  4  )  )  UZX/.A  •l,T,-,=  rfff+,»U£*SeaBfiBB8aaiBBBaaBf«££*£»fiSf^F.feyvAXX  Z  1  1  )  )  4  4  .= 
.  =  =  44  )  )  ll/xy,.A.VcTi=»'4»f£e€IS8aBBBBBaaBSBaB88taa-i‘iil:£8Ht--'£«'..-1A>.>  /  ll  )  )4  1  = 

.  =  ==44 )  1  i/z>  /  a.'  ••A=,!=»  »?i£433*B8aa*aaa»BBBBBEBtaa3-*ai£i4i*Ft--,'-AAxx/  n  1 44  -_  = 

.  =  =  =**)  )i//x)v  it)75=)4i3-aniiaEC8*aB8aBi*a8aB*«**4ic.£'fi<4g=ic -r-'AX// 111)44=- 
.  =  =  =  4  4  )  )  1  1Z  XX  ‘.A  'V-'=Jf4.4?,H4isaeSiBBBBBBBBBBB8SS44S£fK*1F.'-t-  yA  AX  X  i  11)  >  4  4  =.- 
..—-444)11  Z/XXA,'i;',-x,fc;t.4^t4,l**eElBlB8eBBea8SISiiltf£-iJOci-  •  -1AAXZZ  1  )  )4  4  =  -  = 
..=  =  =44)11  1/  XXA  A"-ir.S?,4aii  JiilsaaaaaBfE5iSIf««££f  AAXXZZ1  )  )  4  4=  =  = 

.  .  .  =  =4  4  )  )  11  <./  XX  Af/  v/fr  fef  4  ££££££* |:iaaaaaasc««ijfiiaa'r-b'1  AXXZZ  1  1  )  4  4  4  =  --. 
...=  =  =  4  4  )  )  11  //X AA-  ',ti-F4»'4»  J£-£4434t5***IS49-aia-aS;(Ett^!=ey/AAXXZ  1  1  )  )44  =  =  =  . 
...  =  =  =  +  ♦  4  )  1  IZ  Z  XXAAZl-.^  =  4.)11^ViJ4-jT»£*-»4aaj£J-Jtfffitttfrgf  M=  AAXX  ZZ1  1  )  +  44  =  --  .  . 
,...==  =  *4))ll//XXAAy  ■■=,<-=  f-HA-J  «  V4i  i'tii-U  iftat  A  AXXZZ  11)144  =  =  =.. 

-....  =  =  444 j  j  i  1.-/XXAA-"  .-=='=-ii4  ■•■"  \\yxzzi  1 1)44=  =  =  ... 

-....=  =  =  44)))11  ZZ/XA-*/  -  — - '  V=5=  -4  th+4W  WM  (hl4>4;ij.i#14gfe(--iH  •i"AAXXZZll)))44=  =  =... 
- ....=  ==  44  »  )  )  1  1/./XX  ■*.'>*»  ;t  r)=  =.--JF,-=ft.4(4  4.H>  f  ggg  gf  f  ti  f..-l  A  A  AX  X  Z  Z  1  1  )  )  )  4  4===  - - - 

—  .  .  .  .  =  =  =44  f  )  )  )  1 1  zz  yx>  '.A‘  '=  — r  =  r  =  ^Pi=r=  rittr  VV,‘AAXXX  ZZll)))  +  4  =  =  =....-- 

- .... =  =  =444)1)11/  / /yx/MA-  '4|v  A«.;\  X  X  X  /  Z  /  11)1)444==  =  ....  — 

- ....-=  =  =  4  4  4  4  1  )  1  1  1  1  1  /XX  XX  A"  AAA  ■.  !  •’  '•  7!  AAAAAXXXX/ZZ  l  1  l  )  )  444  4=  =  =  .  .  .  . - 

- .  ...  =  =  =  =  4  4  4  ))  )  )  1  1  !  /  /  7  /  X  a  x  x  A  A  A  A  <’  .XXXXyXZZZ  1  I  1  !  )  1)444  =  -  =  -.... - 

- ...  .=  =  =  =  4444  )))!  11  lZZ//iA<XXXXX/ XXX Z//Z/11U)  1)4444  =  =-  =.,.. - 

- - - -  _  _  4  4  4  4  )  )  )  )  >  1  I)  /  /,  Z  t  UH  ill  Lll  \  1  1  i  )  )  )  )  4  *♦  4=-=  -= - - 

- - '-=  =  -*  4  4  4  4  )))))  1  !  1  1  ;  |  !  1  1  1  1  11  1  1  1  !)))))  44  4  4  4==  =  = - - 

- - ;  -  -  =  4  4  *  4  4  ))))))))  1  111  1  1  1  ))))))))  4  4  4  4  4=  =  =  =  = - - 

- - =  -  =  -  -  -  ***  +  4*  1  |  )  )  )  )  ))  )  )  )  )  ))  )44-4*44  =  ==  =  =  = - - 

- - ..---  =  .-  =  =  4*t**4***444*4**444*4  -  =  =  ==  =...... - 


7  =  =  -  *  4*44*4  *  44  **-  =  : 


SUBMARINE  HEADING:  UNIFORM  DISTRIBUTION  0-27 T 
SUBMARINE  VELOCITY:  1  UNIT/UNIT  OF  TIME 
ELAPSED  TIME:  2.0  UNITS 

NUMBER  OF  SONOBUOYS:  2  NUMBER  OF  "LOOKS" :  5 

TIME  OF  "LOOKS":  1.0,  1.25.  1.5,  1.75,  2.0 

Figure  3-7  PROBABILITY  DENSITY  OF  SUBMARINE'S  POSITION  WITH  KNOWN 
CONSTANT  VELOCITY  AND  UNIFORM  HEADING:  METHOD  IL 
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ELAPSED  TIME:  2.0  UNITS 

NUMBER  OF  SONOBUOYS:  2  NUMBER  OF  "LOOKS":  5 
TIME  OF  'LOOKS":  1.0,  1.25,  1.50,  1.75,  2.0 

Figure  3-8  PROBABILITY  DENSITY  OF  SUBMARINE'S  POSITION  WITH  KNOWN 
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SUBMARINE  HEADING:  UNIFORM  DISTRIBUTION  0-27 T 
SUBMARINE  VELOCITY:  V  =  V  =  N  (0,1) 

ELAPSED  TIME:  2  UNITS  y 

NUMBER  OF  SONOBUOYS:  2  NUMBER  OF  "LOOKS":  5 
TIME  OF  "LOOKS":  1.0,  1.25,  1.5,  1.75,  2.0 

Figure  3-9  PROBABILITY  DENSITY  OF  SUBMARINE'S  POSITION  WITH  RAYLEIGH 
DISTRIBUTION  OF  VELOCITY  AND  UNIFORM  HEADING:  METHOD  I 
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SUBMARINE  HEADING:  UNIFORM  DISTRIBUTION  0-2  71 
SUBMARINE  VELOCITY:  Vx  =  V  =  N  (0,1) 

ELAPSED  TiME:  2  UfviiTS 

NUMBER  OF  SONOBUOYS:  2  NUMBER  OF  "LOOKS":  5 
TIME  OF  "LOOKS":  1.0,  1.25,  1.5,  1.75,  2.0 

Figure  3-10  PROBABILITY  DENSITY  OF  SUBMARINE’S  POSITION  WITH  RAYLEIGH 
DISTRIBUTION  OF  VELOCITY  AND  UNIFORM  HEADING:  METHOD  II 
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SUBMARINE  HEADING:  UNIFORM  DISTRIBUTION  0-2 7T 
SUBMARINE  VELOCITY:  V  -  V  =  N  (0,1) 

ELAPSED  TIME:  2  UNITS  y 

NUMBER  OF  SONOBUOYS:  2  NUMBER  OF  "LOOKS":  5 
TIME  OF  "LOOKS":  1.0,  1.25,  1.5^1.75,  2.0 

Figure  3-11  PROBABILITY  DENSITY  OF  SUBMARINE'S  POSITION  WITH  RAYLEIGH 
DISTRIBUTION  OF  VELOCITY  AND  UNIFORM  HEADING:  METHOD  EL 
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portion  of  the  Search  Director  study  effort  has  been  concentrated  toward  the 
full  development  of  this  computer  program. 

4.  Start  and  End  Grids  Development 

The  initial  probability  density  distribution  of  target  location 
within  a  SPA  at  i-0  is  described  by  joint  normal  probability  distributions. 
The  computer  model  initially  computes  f  [t^)  Ax  and  f  (yy)Ay  and  stores  each 
in  a  one-dimensional  array.  The  model  as  presently  configured  saves  on 
computer  core  size  and  a  minimal  amount  on  computer  run  time  compared  with 
computing  the  above  functio.  s  over  the  total  initial  distribution  (SPA)  and 
then  storing  this  array  in  machine  memory. 

The  locus  of  points  in  the  Z-y  plane  such  that  £(X,  y  )  is  a 
constant  Ci  is  an  ellipse  given  by:  ' 


Therefore,  the  boundary  of  the  starting  grid  in  the  model  is  an  ellipse  which 
is  determined  by  C,  .  The  value  of  is  an  input  to  the  model.  The 
target  submarine  is  thus  assumed  to  be  within  an  ellipse  of  specified  area 
(SPA  SIZE)  with  a  given  confidence,  e.g.,  90%  at  time  =  0.  The  region  outside 
the  ellipse  is  not  considered  in  the  analysis. 

The  Ax  by  A y  start  grid  size  need  not  be  of  equal  size  in 
the  X  and  u  direction.  That  is.  Ax  is  not  required  tc  be  the  same  size 
as  Ay  .  This  allows  f (  X  )  and  f (  y)  to  be  subdivided  into  an  equal  number 
of  grids  even  when  their  variances  are  unequal.  The  end  grid  calculated  at 
the  update  time  t  is  circular  in  shape.  The  grid  center  and  radius  are 
inputs  to  the  model.  A  limitation  of  the  model  is  that  for  the  end  grid  Ax  = 
Ay  ,  the  value  of  AX  should  not  be  larger  than  the  maximum  of  the  values 
of  Ax  or  Atj  of  the  start  grid. 
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5. 


Probability-of-Miss  (PM)  Grid  Development 


One  limitation  of  the  present  computer  model  is  that  all 
deployed  sonobuoys  are  assumed  to  have  the  same  probability  of  detection 
curve.  The  probability  of  detection  of  the  target  by  a  sonobuoy  is  assumed 

_  qc  R  2 

to  be  of  the  form  e  where  R  is  the  distance  from  the  target  to  the 

sonobuoy  and  cc  is  an  input  to  the  model.  All  sonobuoys  being  monitored  are 
assumed  to  be  sampled  at  the  same  time.  It  is  also  ..ssumed  that  the  sampling 
rate  remains  constant  throughout  the  search  time. 

The  probability  that  a  sonobuoy  does  not  detect  the  target  is 

given  by: 


In  the  computer  model  the  PM  grid  is  calculated  utilizing  the  following 
relationship: 

pm{J-  -  1-  exp  /- +(i-  1)  Az)  +  (^j-  +  (j  -  1)  Au) for  C,j  x.  K 

where  c,j  represents  the  number  of  z  grids  and  y  t>rids  from  the  sonobuoy 
to  the  target  and  K  is  an  input  to  the  model.  Also,  the  Ax  A  u  grid  size 
of  the  Pm  grid  can  be  an  integer  subdivision  of  the  end  grid  Ax  Ay  size. 
This  allows  for  a  finer  sampling  of  the  PM  curves.  Hence,  the  model  has 
three  different  grid  sizes  that  can  be  used  and  are  listed  below: 


Start  Grid 

AVsg 

End  Grid 

A*£G 

PM  Grid 

A*  pm 

Also,  <x  may  be  adjusted  to  become  a  function  of  velocity  for  a  limited  number 
of  target  velocity  and  heading  distributions  and  when  the  velocity  is  a  random 
variable.  This  will  adjust  the  probability  density  map  for  these  limited 
cases  to  account  for  an  increase  in  source  level  of  the  target  a:>  it's 
velocity  increases.  The  model  is  presently  being  extended  to  include  any  sensor 
Probability  of  Detection  versus  range  functions  for  both  direct  path  and 
convergence  zone  detection. 
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For  each  "look"  by  Che  sonobuoy  field,  the  model  checks  each 
sonobuoy  to  determine  whether  that  sonobuoy  was  monitored  for  that  "look". 

If  the  sonobuoy  is  not  being  sonitoredj pm,  is  set  equal  to  one.  The 
computer  model  calculates  the  x  and  y  component  of  range  of  each  MC  submarine 
to  each  sonobuoy  for  all  sonobuoys  which  are  being  monitored.  If  either  the 
X  or  y  components  (scaled  to  the  number  of  grids)  is  greater  than  K  ,  -prn-j 
is  again  set  equal  to  one.  If  x  and  ^  are  both  less  than  K  ,  the  model 
looks  up  the  value  of  p m  -  in  the  PM  grid.  The  accumulated  probability  of 
miss  of  the  MC  submarines  is  then  the  product  of  all  the  probability  of  misses 
for  each  sonobuoy  on  each  "look". 

6.  Methods  of  Assigning  Velocity  Components  to  MC  Submarines 

The  objective  in  analyzing  different  methods  of  assigning 
velocity  components  to  MC  submarines  is  to  find  a  method  which  allows  the 
minimur.:  number  of  MC  submarines  to  be  placed  in  each  A#S(.  by  AySQ  grid  in 
order  to  satisfactorily  compute  the  probability  density  map.  Three  methods 
are  available  to  obtain  the  set  of  numbers  representing  the  headings  of  the 
MC  submarines  for  this  distribution: 

I.  Generate  a  set  of  random  numbers  from  this  distribution. 

II.  Equally  space  the  headings  of  the  MC  submarines  between 
0  and  2  fr  and  weight  the  count  of  the  MC  submarines  by 
the  probability  that  the  target  actually  has  that  heading 
within  some  A  0  range. 

III.  Space  the  heading  of  MC  submarines  at  equal  probability 
increments  between  0  and  2  Tf  . 

a.  Uniform  Heading  Case 

The  first  attempt  at  assigning  velocity  components  to  the 
submarines  has  been  to  use  a  random  number  generator  for  the  "random" 
portion  of  the  velocity  and/or  heading  distribution.  An  an  example,  consider 
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the  case  where  the  magnitude  of  the  velocity,  ,  is  known  and  the  heading 
of  the  distribution  is  assumed  to  be  uniformly  distributed  between  0  and  2v 
Let  N  be  the  number  of  MC  submarines  placed  in  each  A  x  by  Ay  grid. 

A  uniform  random  number  generator  supplies  a  list  of  /V  "random"  numbers 
between  0  and  1.  Each  number  of  the  set  is  them  multiplied  by  2  7f  to 
obtain  a  set  of  random  numbers  which  are  uniformly  distributed  between  0  and 
2  7f  .  The  velocity  components  for  the  ith  MC  submarines  are  then  formed 
from  the  following  relationship: 

=  v  cos  (Qc  ) 

Z  i 

7/yt  =  ^  s,n  ) 

where  8 L  is  the  ith  number  in  the  random  number  set. 

Another  approach  is  to  equally  space  the  headings  of  the 
submarines  in  each  grid  around  0  to  2  ff  t  that  is,  t.c  assign  "expected 
values"*  of  headings  to  the  MC  submarines.  Thus,  if  four  MC  submarines  are 
to  be  placed  in  each  Ax  by  Ay  grid,  the  following  headings  would  be 
assigned  to  the  MC  submarines  in  that  grid: 

0,  90”,  180°,  270° 

Sample  result.';  using  these  two  methods  are  presented  in  Figure  5  12. 
Comparison  of  the  two  methods  of  assigning  velocity  components  to  MC 
submarines  is  shewn  in  determining  the  probability  density  distribution  of 
a  moving  target  submarine  having  a  known  constant  velocity  and  uniform 
heading.  In  Curve  I  the  headings  of  the  MC  submarines  are  obtained  utilizing 


*"Expected  value"  is  used  in  the  following  souse:  If  a  probability  density 
function  is  divided  into  hi  subdivisions  in  which  the  area  in  each  segment  is 
equal  (equal  probabilities)  and  N  random  samples  are  obtained  from  the 
distribution,  it  is  expected  that  one  sample  would  fall  in  each  of  the  V  sub 
divisions. 


rii-3i 


UNIFORM  HEADING:  FOR  TWO  METHODS  OF  ASSIGNING  VELOCITY 
COMPONENTS  OF  MC  SUBMARINES 


the  random  number  generator.  In  Curve  II  the  headings  are  equally  spaced 
between  0  and  2  tt  .  Unlike  Curve  I,  Curve  II  corresponds  closely  to  the 
exact  results  presented  in  Reference  1.*  These  results  indicate  that  equally 
spacing  the  headings  of  the  MC  submarines  allows  the  minimum  number  of 
submarines  to  be  used  m  each  grid. 

b.  Pi ased  Heading  Case 

Consider  the  example  where  the  magnitude  of  the  velocity  is 
known  and  the  heading  is  a  random  variable  described  by  the  following 
probability  distribution: 

p(8,  A  m)  '  (a~)  SXP  iAw  C5S  (e)}  ~  V  "  9  '  V 

The  function  lc  (_A.„,  )  is  a  modified  Bessel  function  of  the  first  kind 
which  is  included  so  that  the  density  will  integrate  to  unity.  The  parameter 
Am  can  be  regarded  simply  as  a  parameter  that  controls  the  spread  of  the 
density.  From  the  above  equation  it  is  determined  that  for  Am  =  0. 


As  A.m  increases,  the  density  becomes  more  peaked  and  as  it 

approaches  the  case  where  the  heading  is  known.  Therefore,  the  above 
distribution  allows  a  wide  range  of  headings  distribution  to  be  modeled.  Due 
to  the  kind  of  results  obtained  for  the  uniform  heading  case,  no  attempt 
has  been  made  to  model  a  random  number  generator  for  the  above  distribution. 

Results  for  the  biased  heading  case  utilizing  Methods  II  and 
III  are  presented  in  Figure  3-13.  These  indicate  that  Method  II  (Illustration  13a) 
provides  a  better  representation  of  the  probability  density  map  in  the  fourth 
quadrant  than  Method  III  (Illustration  13b).  However,  the  reverse  situation  occurs  for 


‘Reference  i:  Status  Summary  -  Research  at  CAL  on  Analysis  of  Airborne  Anti- 
Submarine  Warfare  (ASW)  Operations  CAL  Report  -  no  number. 

May  1962,  SECRET. 
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Figure  3-13  PROBABILITY  DENSITY  OF  SUBMARINE'S  POSITION  WITH  CONSTANT 

VELOCITY  AND  RANDOM  HEADING:  FOR  TWO  METHODS  OF  ASSIGNING 
VELOCITY  COMPONENTS  OF  MC  SUBMARINE'S 


the  second  quadrant.  It  is  believed  that  both  methods  can  be  extended  to 
provide  improved  results.  However,  this  has  not  been  attempted  up  to  the 
present  time. 


c.  Random  Velocity  Component  Case 


Consider  the  example  where  the  velocity  is  a  random  variable 
described  by  the  following  joint  normal; 


fUz-  *>;-  z  r 


~v'n 


e*p  1 


The  velocity  components  for  this  distribution  may  be  obtained  from  a  random 
number  generator  in  the  following  way:  The  first  step  is  to  obtain  two  sets 
of  N  random  numbers  from  the  following  >V  (  0,  1  )  distribution: 


_ l_ 

■/It? 


ex  p 


B  *’} 


The  next  step  is  to  multiply  each  number  in  one  of  the  sets  by  <r ^  and  add 

vXq  .  The  third  step  is  co  multiply  each  member  of  the  other  set  by  crv^ 
and  add  Vyo  .  The  final  step  is  to  pair  the  numbers  of  the  two  sets  to 
obtain  a  new  set  which  corresponds  to  a  set  of  random  numbers  from  f  (  v-x,  vy). 
bach  member  of  this  set  is  then  paired  with  one  of  the  MC  submarine  as 
its  velocity  component. 


Another  approach  in  assigning  velocity  components  to  the 
submarines  is  to  obtain  a  set  of  "expected  values"  of  velocity  components 
from  the  distribution  f  (  ,  v-y  )  as  follows: 


Let 


vx  - 


- 


then 


Jo 


f  (V\ 


Ufcr^cj^ 


1 11-35 


Let 


v  =  v.  - 

Y  9  cr„ 


then 


Wow  let 


'<*-*>’  17^  exp\- ^ 

*x  L  J 


=  7^ 


v  =  i  TS-*  +  'v-  3 


The  random  variable  ir  has  a  Rayleigh  density  given  by: 

2r  \  2r 2  1 

~  e*P\~  7^— 1 

rt 


3  =  -jrs-  e*p 


2<V,  j 


for  zr  >  0 


Solve  for  M  velocities  along  the  curve  g  (  V'  )  spaced  at  equal  probability 
increments. 


Next,  take  the  first  of  the  A?  values  of  vc-locity,  z/^  ,  above 

and  divide  0  into  equal  increments  over  2  V  in  the  same  manner  as  in  the 
case  for  uniform  heading.  Tne  reason  for_this  is  that  the  random  variable 

a ,  Un-’ 

is  uniformly  distributed  between  0  and  2  tt  . 

v:z(6l)=  IP  cos(6{) 

!Py  (3-)  -  IP  j/n  ( Qc  ) 

It  is  now  necessary  to  transform  the  set  of  expected  values  from  the 
distribution  f  ( ir  ,  V-y  )  to  a  set  from  /  ( ,  vu  ). 

Therefore, 

*i  <s‘  >  •  Tp  “y  (e‘>  *  \ 

if 

v-x(dL)  *  vx(dc)  +  vXq 


in-36 


Let  ex:  be  the  angle  between  of  the  vy  vy  plane  and  x,  of  the  Zy  SPA. 

Then  rotate  the  values  of  vx  (  6;  )  and  vy  (  )  into  the  % y  plane  of  the  SPA 

by  the  following  equations: 

V*(%  )  -  cos(ct)  +  vy(6L)  cos  (<*  +  jr) 

=  v-z(9l)  sm(oc)  +  S/n  («  +  4r) 

The  above  operation  is  repeated  for  each  of  the  M  original  velocities  to 
obtain  the  complete  set  of  velocity  components  to  model  f  (  vy  ,  vy  ) . 

Sample  results  are  presented  in  Figure  3-14  op  the  probability 
density  of  a  submarine's  position  for  a  Rayleigh  distribution  in  velocity  and 
uniform  heading  utilizing  the  random  number  generator  and  expected  value 
methods.  This  is  also  presented  for  the  exact  solution  for  comparison 
purposes.  In  the  three  illustrations  (14a,  14b  and  l4c)  and  tfn,  are  both 
described  by  a  normal  distribution  with  a  mean  of  zero  and  a  standard 
deviation  of  12  knots.  The  resultant  vector  velocity  is  described  by  a 
Rayleigh  density  and  hence  the  expected  value  of  the  vector  velocity  is 
approximately  15  knots.  The  initial  SPA  has  a  standard  deviation  of 
25  n.  miles  for  both  the  z  and  y  components.  The  A*  by  A  y  grid  size  is 
4  x  4  n.  miles.  The  display  time  is  8  hours. 

Illustration  14a  is  the  exact  solution  density  map  for  the 
sample  problem.  In  illustration  14b,  the  velocity  components  are  generated 
from  a  normal  random  number  generator.  The  velocity  components  for 
illustration  14c  have  been  obtained  using  expected  values.  The  number  of  MC 
submarines  in  each  grid  are  the  same  for  both  illustrations  14b  and  14c. 

These  examples  indicate  that  using  the  expected  value  method  gives  the  best 
results  since  they  more  closely  approximate  those  shewn  for  the  exact 
solution. 


Additional  results  are  presented  in  Figure  3-15  using  the 
expected  value  method.  The  parameter  values  used  are  the  same  as  those  for 
Figure  3-14  except  for  the  number  of  MC  submarines  placed  in  each  grid. 
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Figure  3-14  PROBABILITY  DENSITY  OF  SUBMARINE'S  POSITION  WITH  RAYLEIGH 
DISTRIBUTION  IN  VELOCITY  AND  UNIFORM  HEADING:  FOR  TWO 
METHODS  OF  ASSIGNING  VELOCITY  COMPONENTS  TO  MC  SUBMARINES 


Illustration  15a  in  Figure  3-15  is  the  result  of  a  computer  run  in  which  the 
function  a  (  £r  )  is  sampled  at  a  relatively  small  sample  size,  M  ,  and  a 
relatively  large  number  of  MC  submarines  is  assigned  to  each  sampled  velocity. 
In  15b.  g  (  v  )  is  sampled  at  a  large  sample  size,  M  ,  but  the  number  of 
submarines  assigned  to  each  sampled  velocity  is  small. 

Referring  to  15a,  a  readily  discerned  effect  resulting  from 
under  sampling  cy  (  W)  is  that  a  "hole"  in  the  probability  density  pattern 
appears  at  the  datum.  On  the  other  hand,  assigning  relatively  few 
submarines  to  each  sampled  velocity  (15b)  causes  the  constant  contour  levels 
to  have  the  effect  of  a  circle  modulated  by  a  sine  function.  Thus,  additional 
tradeoff  analyses  are  necessary  to  minimize  or  eliminate  these  adverse 
effects. 


7.  Sub-Optimal  Sonobuoy  Drop  Point  Determination 

A  limited  effort  has  been  conducted  on  a  modification  to  the 
Monte  Carlo  program  for  use  as  a  possible  improved  aid  in  determining  sub- 
optimal  sonobuoy  drop  points.  Instead  of  solving  for  the  joint  probability 
that  the  target  is  in  some  grid  at  time  T  and  that  the  target  will  not  be 
detected  in  moving  to  this  grid,  the  modified  program  calculates  the  joint 
probability  that  the  target  passes  through  some  grid  during  the  search  time 
from  to  Tz  and  that  the  target  is  not  detected  during  this  time  by  the 
deployed  sonobuoys. 

The  modified  program  would  be  used  in  the  following  manner: 
First,  the  operator  determines  the  time  on  s'ation  (from  7"f  to  Tz  ) .  The 
model  then  computes  the  probability  that  the  target  will  pass  through  each 
A  %  by  Ay  grid  during  this  time  period.  Having  this  display  the 
operator  then  decides  on  a  drop  point  for  one  or  more  sonobuoys.  The  model 
then  calculates  the  joint  probability  that  the  target  passes  through  the 
A  x  by  Ay  grid  and  is  not  detected  by  the  deployed  sonobuoys.  For  this 
calculation  the  operator  must  make  a  reasonable  estimate  of  the  monitoring 
schedule  of  the  deployed  sonobuoys. 
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Presently,  a  problem  area  exists  in  this  model  which  is 
caused  by  the  use  of  a  rectangular  coordinate  system.  The  reason  is  that 
submarines  with  a  heading  of  0°  may  be  counted  as  passing  through  more  grids 
than  those  with  a  heading  of  45°  for  example.  A  correction  factor  has  been 
inserted  into  the  model  to  account  for  this  difference.  However,  a  complete 
analytical  justification  of  this  factor  has  not  vet  been  made. 

Sample  results  from  this  model  are  presented  in  Figure  3-16. 
Illustrations  16a  and  16b  respectively  show  the  probability  density  maps 
without  and  with  application  of  the  correction  factor.  In  both  cases  the 
target’s  velocity  of  13  knots  with  a  uniform  heading  between  0  and  2  w  is 
assumed  to  be  known.  The  grid  size  is  2.5  x  2.5  n.  mi.  The  SPA  is  described 
by  a  joint  circular  normal  distribution  with  a  standard  deviation  of 
25  n.  mi.  In  summary,  although  complete  analytical  justification  for  the 
correction  factor  has  not  been  made  computer  run  results  (16b)  indicate  that 
it  appears  to  work  satisfactorily. 

8.  Analog-Hybrid  Computer  Program 

A  limited  effort  has  also  been  conducted  toward  determining  the 
feasibility  of  using  an  analog/hybrid  computer  system  to  obtain  the  probability 
density  maps  of  target  uncertainty  areas.  This  limited  program  has  been 
completed  on  the  CAL  COMCOR  Cl  5000  Analog  computer.  Typical  results  of  the 
CRT  display  obtained  from  this  program  are  illustrated  in  Figure  3-17.  The 
tactical  situation  shown  is  for  the  case  where  the  velocity  of  the  submarine 
is  a  constant  value  and  it's  heading  is  uniformly  distributed  between  0  and 
2  rr  .  Referring  to  Illustration  17a,  one  "look"  is  taken  with  two  sonobuoys 
offset  on  both  sides  of  the  center  of  the  uncertainty  area.  In  17b,  two 
looks  with  one  sonobuoy  are  shown  for  the  case  in  which  the  sonobuoy  is 
dropped  at  the  center  of  the  uncertainty  area. 

A  general  purpose  computer  could  be  used  in  conjunction  with 
the  Analog  computer  (Hybrid  system).  This  system  would  use  the  GP  computer 
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Figure  3-16  ACCUMULATED  PROBABILITY  DENSITY  OF  SUBMARINE'S  POSITION  WITH 
KNOWN  CONSTANT  VELOCITY  AND  UNIFORM  HEADING:  EFFECT  OF 
CORRECTION  FACTOR 
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primarily  for  storage  while  the  analog  computer  would  be  used  to  "compute" 
the  probability  density  maps. 

9,  Possible  Search  Director  Display  Techniques 

A  limited  study  has  been  conducted  of  several  techniques  for 
displaying  probability  density  maps  of  target  uncertainty  areas  in  which 
sonobuoys  have  been  deployed.  Figure  3-18  illustrates  several  display 
technique  which  have  been  photographed  from  a  CRT  connected  to  the  output 
processor  of  the  CAL  analog  computer.  These  presentations  vary  the  amount  of 
fill-in  between  the  calculated  value  for  each  Az  by  Ay  grid.  In 
i’ lustration  18a  the  fill-in  is  almost  complete.  In  18b  a  partial  fill-in 
is  utilized  for  display  purposes.  In  18c  the  probability  level  is  indicated 
both  by  the  intensity  of  the  point  and  by  its  size. 

A  capability  also  exists  in  future  study,  for  the  analog 
computer  to  pass  the  output  grid  through  a  non- linear  amplifier.  This 
capability  provides  a  method  for  investigating  the  effects,  for  example,  of 
intensifying  the  highest  probability  density  levels, 

10.  Comparison  of  Display  Techniques  Utilizing  the  CAL  Flying 

Spot  Scanner 

A  limited  effort  has  been  conducted  on  the  CAL  Flying  Spot 
Scanner  for  the  purpose  of  comparing  various  display  techniques.  A  description 
of  this  CAL  equipment  is  summarized  in  Appendix  A.  In  brief,  the  values  cf 
the  calculated  grid  of  the  Monte  Carlo  program  are  converted  to  integer  values 
ranging  from  0  to  60.  The  lowest  number  in  the  grid  is  set  equal  to  zero, 
the  highest  to  60  and  the  intermittent  values  are  linearly  spaced  between  0 
and  60.  A  10/4  x  1024  grid  is  then  formed  from  the  calculated  grid  by 
linearly  smoothing  between  the  calculated  points.  A  color  code  is  assigned  to 
each  of  the  61  levels  in  the  grid.  The  possible  shades  of  gray  for  a  black 
and  white  photograph  are  shown  in  Figure  3-19. 
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LIGHTER  AREAS  INDICATE  HIGHER  PROBABILITY  DENSITY  AREAS 
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Figures  3-20  and  3-21  are  examples  in  which  all  the  different 
levels  of  gray  have  been  used.  The  brighter  areas  indicate  higher  probability 
levels.  Figure  3-20  is  presented  for  the  case  where  the  submarine  target 
velocity  is  known  to  be  12  knots  and  it's  heading  is  uniformly  distributed 
between  0  and  2  ir  .  The  SPA  is  described  by  a  circular  normal  distribution 
with  a  standard  deviation  of  25  n.  mi.  The  displayed  circle  has  a  diameter 
of  252  n.  mi.  The  display  time  for  this  example  is  9  hours.  A  total  of  15 
sonobuoys  are  deployed.  Hach  sonobuoy  is  monitored  on  every  other  "look" 
starting  at  a  time  late  equal  to  5  hours  in  increments  of  15  minutes  for  a 
total  of  24  "looks".  This  photograph  is  for  the  same  probability  density 
map  as  shown  in  Figure  4-7c  of  Section  IV  of  this  report  in  which  an  on-line 
computer  nrintei  lias  been  used  to  create  the  map.  Figure  4-7  also  shows  the 
locations  of  the  deployed  sonGbuovs. 

An  example  is  shown  in  Figure  3-21  in  which  the  velocity  of  the 
target  is  described  by  a  Rayleigh  distribution  with  a  mean  velocity  of 
15  knots.  The  diameter  of  the  display  is  420  n.  mi.  and  the  center  of  the 
figure  is  the  datum  point.  Figure  4- 1 2b  is  a  presentation  of  the  same  display 
using  the  on-line  printer  (Refer  to  Section  IV). 

Figures  5-22  and  3-23  present  the  same  probability  density  map 
that  is  shown  in  Figure  3-20  except  for  the  coding  of  the  shades  of  gray.  In 
Figure  3-22,  31  intensity  levels  have  been  used  which  are  equally  spaced 
over  the  gray-scale  range  of  0-60.  Figure  3-23  shows  six  intensity  levels 
linearly  spaced  over  the  dynamic  range  0-60. 

Figure  3-24  presents  the  same  probability  density  map  that  is 
shown  in  Figure  3-20  except  that  in  this  example,  the  blacker  areas  indicate 
higher  probability  levels.  This  is  accomplished  by  making  a  positive  of 
Figure  3-20. 
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Figure  3-22  FLYING  SPOT  SCANNER  DISPLAY  OF  PROBABILITY  DENSITY  MAP: 
EXAMPLE  3 
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Figure  3-23  FLYING  SPOT  SCANNER  DISPLAY  OF  PROBABILITY  DENSITY  MAP: 
EXAMPLE  4 
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Figure  5-25  shows  an  example  of  a  color  display  with  Figure 
5-26  being  its  equivalent  black-and-white  photograph  using  61  shades  of  gray. 

It  should  be  noted  that  no  effort  has  been  attempted  to  determine  the  "best" 
color  code  and,  therefore,  Figure  5-25  is  presented  only  as  an  example  of  the 
capability  of  the  CAL  Flying  Spot  Scanner. 

Another  method  of  presenting  a  probability  density  map  to  an 
ASK  operator  is  to  display  it  as  a  series  of  contour  lines  representing  iso¬ 
probability  density  levels.  Thus,  in  order  to  evaluate  the  merits  of  this 
method  of  display,  a  contour  plotting  routine  has  been  added  to  the  Monte 
Carlo  model.  Sample  contour  riots  based  on  computer  runs  made  utilizing  this 
routine  are  shown  in  Figures  5-27,  5-28  and  5-29. 

A  contour  plot  of  probability  density  levels  for  the  initial 
SPA  ( N(0,i ) )  is  shown  in  Figure  5-27  at  zero  elapsed  time.  In  Figures  5-28 
and  5-29,  the  target  is  assumed  to  have  a  known  constant  velocity  and  the 
heading  is  uniformly  distributed  between  0  and  2  7T  .  Contour  plots  for 
display  elapsed  times  of  1  unit  and  2  units  are  presented  respectively  in 
Figures  5-28  and  5-29.  The  probability  density  levels  associated  with  zero 
and  one  unit  of  elapsed  time  are  listed  in  Table  5-1.  Those  for  two  units  of 
elapsed  time  are  listed  in  Table  5-2. 

In  summary,  the  limited  results  achieved  to  date  provide  some 
indication  that  an  intensity  display  in  shades  of  gray  can  supply  ASK  aircraft 
personnel  with  the  necessary  information  required  for  conducting  Search 
Director  operations.  Still  to  be  finally  determined,  however,  are  the  number 
of  intensity  levels  needed  and  the  size  of  the  display  scope. 
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TABLE  3-1  PROBABILITY  DENSITY  LEVELS  FOR  ZERO  AND  ONE  UNIT  OF  ELAPSED  TIME 


PROBABILITY  DENSITY  LEVEL 

li(u,v)  u 

V 

0 

0  ^  Prob . 

<  .0005 

1 

.0005 

.0001 

2 

.0001 

.0002 

3 

.0002 
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4 
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.0004 

5 

.0004 

.0005 

6 

.0005 

.0006 

7 

.0006 

.0007 
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.0007 

.0008 
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.0008 

.0009 

10 

.0009 

.0010 

11 

.0010 

.0011 

12 

.0011 

.0012 

13 

.0012 

.0013 

14 

.0013 

.0014 

15 

.0014 

.0015 

16 

.0015 

:  3-2  PROBABILITY  DENSITY  LEVELS 

PROBABILITY  DENSITY  LEVEL 

FOR  TWO  UNITS 

hfu, v) 

OF  ELAPSED  TIME 

u  V 
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0  £ 

Prob.  ^  .0002 
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.00028 
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GENERAL  STUDY  RESULTS 


1 .  Comparison  of  Sample  Results  for  Exact  Solution  and  Monte 
Carlo  Program 


An  illustrative  comparison  of  results  obtained  from  an  "exact" 
solution  versus  those  from  the  Monte  Carlo  program  is  shown  in  Figures  4-1 
and  4-2.  In  these  examples  it  is  assumed  that  the  target's  velocity  components 
are  described  by  //(  o,  1  ) .  The  SPA  is  assumed  to  be  described  by 
a  joint  circular  normal  '.-5th  unity  variance  and  zero  mean.  The  display  time 
is  one  unit.  One  sonobuoy  is  deployed  in  the  uncertainty  area  and  one  "look" 
is  obtained  at  0.75  units  elapsed  time. 

The  similarity  between  the  piobability  density  maps  shown  in 
Figures  4-1  and  4-2  appears  to  be  excellent.  It  should  be  noted  that  the 
computing  time  on  the  IBM  370-165  computer  utilizing  the  exact  solution  is, 
although  not  optimum,  92  seconds  compared  to  less  than  2.5  seconds  for  the 
Monte  Carlo  program. 

2.  Effect  of  Varying  Start  Gtid  Size 

The  results  shown  in  Figures  4-3  and  4-4  indicate  the  effect 
of  varying  the  size  of  the  start  grid.  That  is,  the  effect  of  varying 
to  determine  the  boundary  of  the  start  grid  from  the  following  equation: 


The  importance  of  limiting  the  size  of  the  start  grid  is  that 
doubling  the  size  of  the  grid  (with  Ax  and  Ay  remaining  constant) 
increases  the  computer  run  time  of  the  Monte  Carlo  program  by  more  than  a 
factor  of  four. 
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SUBMARINE  HEADING:  UNIFORM  DISTRIBUTION  0-2  7/ 
SUBMARINE  VELOCITY:  x  &  y  COMPONENTS  NORMAL  WITH 

ZERO  MEAN  AND  UNITY  VARIANCE 


ELAPSED  TIME: 
NUMBER  OF 
S0N0BU0YS: 

TIME  OF  "LOOK": 


1  UNIT 

1,  NUMBER  OF  "LOOKS":  1 
0.75 


Figure 4-1  PROBABILITY  DENSITY  OF  SUBMARINE'S  POSITION  WITH  RAYLEIGH 
DISTRIBUTION  IN  VELOCITY  AND  UNIFORM  HEADING:  EXACT 
COMPUTATIONAL  METHOD 
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SUBMARINE  HEADING: 
SUBMARINE  VELOCITY: 

ELAPSED  TIME: 

NUMBER  OF 
SONOBUOYS: 

TIME  OF  "LOOK'  : 


UN’rORM  DISTRIBUTION  0-2  V 
a  &  y  COMPONENTS  NORMAL  WITH 
ZERO  MEAN  AND  UNITY  VARIANCE 
1  UNIT 

1 ,  NUMBER  OF  "LOOKS":  1 
0.75 


Figure 4-2  PROBABILITY  DENSITY  OF  SUBMARINE'S  POSITION  WITH  RAYLEIGH 
DISTRIBUTION  IN  VELOCITY  AND  UNIFORM  HEADING:  MONTE 
CARLO  METHOD 
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IV-  4 


Figure 4-3  PROBABILITY  DENSITY  OF  SUBMARINE'S  POSITION  WITH  KNOWN 
CONSTANT  VELOCITY  AND  UNIFORM  HEADING:  FOR  VARIATIONS 
IN  START  GRID  SIZE 


Figure  4-4  PROBABILITY  DENSITY  OF  SUBMARINE  POSITION  WITH  RAYLEIGH 

DISTRIBUTION  IN  VELOCITY  AND  UNIFORM  HEADING:  FOR  VARIATION 
!N  START  GRID  SIZE 


For  both  Figures  4-5  and  4-4,  the  SPA  is  assmsed  to  be  a 
joint  circular  normal  with  a  standard  deviation  of  25  n.  mi.  In  Figure  4-5 
the  target  is  assumed  to  have  a  known  constant  velocity  of  12.5  knots  and  a 
uniformly  distributed  herding  between  0  and  2  tr  .  In  Figure  4-4  the  target's 
velocity  is  described  by  a  Rayleigh  distribution  having  a  mean  of  15  knots 
and  a  uni  fora  heading  between  0  and  2 n  . 

The  following  list  associates  the  radius,  C?  ,  of  the  start 
grids  with  each  of  the  illustrations  shown  on  the  two  Figures  4-5  and  4-4. 

1 1  lustration  £, 


4- 5a 

1  cr 

4 -5b 

2  <r 

4-5c 

2.4  cr 

4-4a 

1.6c- 

4-4b 

2<r 

4-4c 

5<r 

From  these  limited  results,  the  2  cr  boundary  appears  to  be  satisfactory  for 
use  in  analyses.  Thus,  for  the  additional  results  presented  in  this  report, 
the  2cr  boundary  is  used  for  the  start  grid.  This  encompasses  approximately 
90%  of  the  probability  density  which  describes  the  SPA. 

5.  Effect  of  Varying  the  Variance  in  Target  Heading 

Figure  4-5  presents  results  obtained  due  to  varying  the  amount 
of  spread  (variance)  of  the  probability  density  function  describing  the 
target's  heading.  In  this  example  the  target's  velocity  is  given  by  the 
following  probability  density: 

j"  % 

f(v)  =K  exp  j  ^  |  -1b  v-  <  16 
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In  this  case  the  velocity  is  assumed  to  be  a  constant  (12.5 
knots)  plus  a  random  part  described  by  a  normal  distribution  over  the 
interval  from  -2 cr  to  2  cr  (where  cr2  -  64  ).  The  constant,  K,  is  inserted 

so  that  the  function  will  integrate  to  one.  It  should  be  noted  that  only  the 
nagnitude  of  the  velocity  is  determined  from  ~  (  y  ) . 

Tne  heading  distribution  is  given  by: 

p  <s-  s  Tb-TTZJ  exp  {A"  "s  5 '} 

Tne  results  shown  in  Figure  4-5  are  obtained  for  a  number  of  values  of  x\m  . 
Tne  values  of  A.m  associated  with  each  of  the  illustrations  shown  are  listed 
below: 


Illustration  Am 


4-5a  1 

5b  2 

5c  4 

Sd  8 


Curves  of  p  (  6  ,  Am  )  can  be  found  in  Reference  2*  page  538 
for  a  number  of  values  of  A  m  . 


4.  Effect  of  Varying  the  Variance  in  Target  Velocity 


A  presentation  of  selected  examples  to  show  the  effect  of 
varying  the  variance  of  target  velocity  when  the  target's  heading  is 
uniformly  distributed  between  0  and  2  is  presented  in  Figure  4-6. 


‘Reference  2.  Detection,  Estimation  and  Modulation  Theory,  Part  I,  Harry  L. 
Van  Trees,  Published  1968  by  John  Wiley  and  Sons 
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VELOCITY:  k  exp  (-1/2 


SPA  SIZE: 


HEADING: 


(-V-  125  ^2)  -2a<v<20 

(T  2  v  v 

UNIFORM  DISTRIBUTION  0-27T 


°x  =  25  n.mi. 
°y  =  25  n.mi. 
Ax  =  Ay  =  4  n.mi. 


DISPLAY  TIME:  8  HOURS 


Figure 4-6  PROBABILITY  DENSITY  OF  SUBMARINE'S  POSITION  WITH  RANDOM 
VELOCITY  AND  UNIFORM  HEADING:  FOR  VARIANCE  VARIATION  IN 
VELOCITY  DISTRIBUTION 


Tne  probability  density  which  describes  the  velocity  is: 


■p  (v)  - 


l.V-  12.5) 


2  < 


~2<?-r  V  <■  ^°tr 


Tne  target's  heading  is  assured  to  be  uniformly  distributed  between  0  and  2  «*  . 
The  values  of  <rj  associated  with  eadi  of  the  illustrations  shown  in 
Figure  4-6  are  listed  below: 


1 1  lustration 


4-6a  1 

6b  4 

6c  16 

6d  56 


it  is  interesting  to  note  that  as  cr ■  x  increases  the  high 

probability  areas  spread  more  toward  the  datum  than  away  from  it.  That  is, 

the  radius  of  a  circle  contouring  the  outer  edge  of  the  darker  area  has  very 

little  change  over  the  range  of  <rv.z  while  the  radius  of  a  circle  containing 

the  inner  edee  of  the  darker  area  decreases  as  cr *  increases. 

* 

5.  Effect  of  Varying  the  Number  of  MC  Submarines/Start  Grid 


Sample  computer  results  obtained  by  varying  the  number  of  MC 
submarines  in  each  Ax  by  Ay  area  of  the  start  grid  are  shown  in  Figure  4-7. 
This  example  is  for  the  case  of  known  target  velocity  with  a  uniform  heading 
over  the  interval  0-2 if  .  There  are  15  sonobuoys  deployed  in  this  example. 

Each  sonobucy  is  monitored  on  every  other  "look"  period.  The  sonobuoys 
"looks"  are  spaced  at  15  minute  intervals  over  the  period  from  3  to  9  hours 
giving  a  total  cf  24  "looks". 

The  importance  of  limiting  the  number  of  MC  submarines/grid 
is  that  the  computing  time  on  the  IBM  370-165  computer  increases  linearly  as 
the  number  of  MC  submarines  in  each  Ax  by  Ay  grid  increases.  The  number 
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NUMBER  OF  MC 
SUBMARINES/GRID  =  12 
7a 


NUMBER  OF  MC 
SUBMARINES/GRID  =  18 
7b 


of  '1C  submarines/grid  associated  with  each  of  the  illustrations  shown  in 
Figure  4-7  is  summarized  below: 

Illustration  MC  Subparines/Grid 

4- 7a  12 

7o  18 

7c  24 

7d  40 

In  comparing  Illustration  4-7a  (12  .MC  submarines/grid)  with  4-7d  (40  MC 
submarines/Grid)  it  is  noted  that  the  probability  densities  resulting  from 
the  number  of  MC  submarines/grid  are  dependent  on  the  location  of  the 
sonobuoys.  It  is  for  the  purpose  of  showing  this  effect  that  the  sonobuoys 
are  placed  in  an  offset  pattern  with  respect  to  datum. 

The  computing  time  on  the  IBM  370-165  computer  (central 
processing  unit  time)  for  Illustration  4-7d  (40  MC  subraarines/grid)  is  4.1 
seconds  while  the  computing  time  for  4-7c  (24  MC  submarines/grid)  is  2.7 
seconds.  These  times  are  for  the  case  in  which  the  size  of  the  start  grid  is 
determined  by  a  2  cr  boundary.  If  the  range  of  the  target  to  the  sonobuoy 
is  greater  than  12  n.  mi.,  the  probability  that  the  sonobuoy  detects  the 
target  is  assumed  to  be  zero.  An  increase  in  the  values  of  either  of  the 
above  parameters  will  cause  an  increase  in  the  computing  time. 

6.  Iiffect  of  Varying  Search  Aircraft  Time  Late  in  Arriving  at  Datum 


The  effect  of  varying  ASW  search  aircraft  time  late  in 
arriving  at  datum  of  a  SPA  area  and  conducting  searui  operations  for  a 
constant  6  hour  period  is  presented  in  Figure  4-8.  In  this  example,  the 
velocity  of  the  target  is  described  by  the  following  joint  normal  probability 


distribution: 

f  (*x  >  *s} 


TV- 12 


DISPLAY  TIME:  11  HOURS 
TIME  LATE:  5  HOURS 

8c 

O  DATUM 

2  2 

VELOCITY  ~  — —  exp  \  -1/2  (-*—  +  [ 

7T  -  96  «  144  16  ' 


NUMBER  OF  SONOBUOYS: 
NUMBER  OF  LUOKS: 


DISPLAY  TIME:  13  HOURS 


TIME  LATE: 

8d 

7  HOURS 

SPA  SIZE: 

°x 

ffy 

=  25  n.mi. 

=  25  n.mi. 

Ax-  Ay  =  5  n.mi. 
15 
24 


Figure  4-8  PROBABILITY  DENSITY  OF  SUBMARINE'S  POSITION  WITH  X  AND  Y 

VELOCITY  COMPONENTS  NORMALLY  DISTRIBUTED:  EFFECT  OF  TIME 
LATE 
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Also,  there  is  an  angle  of  45°  measured  in  the  clockwise  direction  between 
the  i.  -  (j  axis  of  the  SPA  plane  and  the  ^-z'y  axis  of  the  velocity  plane. 
The  SPA  is  described  by  a  joint  circular  normal  distribution  with  a  standard 
deviation  of  25  n.  mi. 

Illustration  4-8a  shows  the  "growth"  of  the  SPA,  without  any 
sonobuoys  deployed  for  an  elapsed  time  interval  of  11  hours.  Time  late  for 
the  other  illustrations  shown  in  Figure  4-8  are  listed  below: 

1 1  lustration  Time  Late,  Hours 

4-8b  3 

4-8c  5 

4-8d  7 


The  results  indicate  the  varying  probability  densities  that  occur  as  time 
late  is  increased  from  3  and  5  hours  tc  a  maximum  of  7  hours. 


7.  Effect  of  Varying  AS1V  Aircraft  Search  Time 


Computer  results  which  show  the  effect  of  varying  ASW  aircraft 
search  time  for  a  selected  time  late  to  datum  are  presented  in  Figure  4-9. 

The  SPA  described  by  a  joint  circular  normal  with  a  standard  deviation  of 
25  n.  mi.  The  velocity  components  of  the  target  are  given  by  the  following 
distribution: 

f  ,  Vy)  -  exP 


Also,  there  is  an  angle  of  45°  measured  from  the  SPA  axis  in  a  clockwise 
direction  between  the  x-  y  axis  of  the  SPA  and  the  axis  of  the 
target  velocity  plane. 


For  each  case  shown  in  Figure  4-9  the  search  operation  begins 
at  seven  hours  time  late.  The  display  times  and  number  of  sonobuoy  "looks" 
associated  with  the  illustrations  as  shown  in  Figure  4-9  are  listed  below: 
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Figure  4-9  PROBABILITY  DENSITY  OF  SUBMARINE'S  POSITION  WITH  X  AND  Y 

VELOCITY  COMPONENTS  NORMALLY  DISTRIBUTED:  FOR  SELECTED 
SEARCH  TIMES 


Illustration 


Display  Tine.  Hours 

-  .  A  "...  * 


Xunber  of  Sonobuoy  Looks 


4 -9a  9  S 

4-9b  11  16 

4-9c  13  24 

Th  results  indicate  that  the  sonobuoys  in  effect  search  out  an  increasing 
portion  of  the  SPA  as  the  search  tine  and  number  of  ’'looks"  by  the  sonobuoy 
field  increase  from  9  to  15  hours  and  froa  8  to  24  "looks"  respectively. 

S.  Effect  of  Linear  Smoothing  Between  Grid  Points 


Figures  4-10  and  4-11  indicate  the  effect  cf  linear  smoothing 
between  calculated  grid  points.  Figure  4-10  presents  results  for  a  medium 
speed  (12  knots)  target  and  Figure  4-11  for  a  nigh  speed  (24  knots)  target. 

The  target's  heading  m  both  cases  is  uniformly  distributed  between  0  and  Zxr  . 
The  SPA  is  described  by  a  joint  circular  normal  with  ?.  standard  deviation  of 
25  n.  mi.  A  field  of  15  sonobuoys  is  deployed  within  the  SPA.  Each 
sonobuoy  is  monitored  on  every  ether  "look".  The  "looks"  are  taken  during 
the  elapsed  time  interval  from  3  to  9  hours  at  increments  of  15  minutes. 

In  Illustration  4- 10b  the  grid  is  calculated  for  a  Ax  =  Ay  = 
6  n.  mi.  Additional  points  are  then  calculated  between  these  grid  points  by 
linearly  averaging  between  adjacent  grid  points.  This  provides  an  effective 
grid  size  of  3  n.  mi.  on  the  display.  Illustration  4- 10a  is  calculated  for 
a  Ax  =  Ay  -  3  n.  mi.  for  comparison  purposes. 

Referring  to  Illustration  4- Jib  in  Figure  11,  the  display  is 
first  calculated  for  a  grid  size  [Ax  by  Ay  )  of  12  n.  mi.  and  then 
reduced  to  one  of  6  n.  mi.  by  the  above  smoothing  technique.  Illustration 
4- 1 1  a  is  calculated  for  Ax  =  Ay  =  6  n.  mi.  in  order  to  cor*  are  the  two 
cases.  The  resulting  probability  densities  for  these  cases  are  fairly 
similar.  Thus,  the  importance  of  the  linear  smoothing  technique  is  in 
reducing  computing  time  by  a  factor  of  approximately  4  for  a  specified  Ay 
by  Ax  grid  selected  for  the  display. 
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Figure  4-10  PROBABILITY  DENSITY  OF  SUBMARINE'S  POSITION  WITH  KNOWN  MEDIUM 
VELOCITY  AND  UNIFORM  HEADING:  EFFECT  OF  LINEAR  SMOOTHING 
BETWEEN  CALCULATED  GRID  POINTS 
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Figure  4-11  PROBABILITY  DENSITY  OF  SUBMARINE'S  POSITION  WITH  KNOWN  HIGH 
VELOCITY  AND  UNIFORM  HEADING:  EFFECT  OF  LINEAR  SMOOTHING 
BETWEEN  CALCULATED  GRID  POINTS 
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9.  Effrct  of  a  Large  Variation  in  Number  of  MC  Submarines/Grid  for 
a  Uniform  Target  Heading  and  Rayleigh  Distribution  in  Velocity 

Figure  4-12  presents  computer  results  obtained  when  target 
velocity  is  described  by  a  Rayleigh  distribution  with  a  mean  of  15  knots  and 
target  heading  is  described  by  a  uniform  distribution  over  the  interval  0-27 7  . 
Illustration  4- 12a  in  Figure  4-12  presents  results  obtained  for  a  display  time 
of  9  hours  for  the  case  in  which  no  sonobuoys  are  deployed.  Illustrations  4-12b 
and  4- 12c  present  results  for  the  case  in  which  15  sonobuoys  are  deployed. 
Sonobuoy  "looks"  are  taken  at  15  minute  intervals  over  the  elapsed  time 
interval  from  3  to  9  hours  for  a  total  of  24  "looks".  Each  sonobuoy  is 
monitored  on  every  other  "look".  The  number  of  MC  submarines /grid  is 
relatively  small  for  the  case  shown  in  Illustration  4-12b  whereas  it  is  large 
for  the  case  shown  in  Illustration  4- 12c. 

The  computing  time  on  the  IBM  370-165  computer  for  the  case 
shown  in  Illustration  4- 12b  is  6.53  seconds  whereas,  it  is  11.64  seconds  for 
4- 12c.  The  boundary  of  the  start  grid  has  a  radius  of  2  cr  .  If  the  target 
is  more  than  10  n.  mi.  from  a  sonobuoy,  the  probability  of  detection  is 
assumed  to  be  zero.  Increasing  the  sonobuoy  detection  range  from  10  n.  mi. 
to  25  n.  mi.  increases  the  computing  time  of  4-12c  from  11.64  seconds  to 
14.75  seconds. 


1 0 .  Effect  of  Changes  in  Target  Velocity  on  Probability  of  Detection 

Figure  4-13  presents  computer  results  for  the  case  in  which 
the  probability  of  detection  of  the  target  by  the  sonobuoy  field  is  a  function 
of  the  target's  velocity.  As  stated  previously  in  Section  III  of  this  Repoit, 
the  probability  of  detection  is  assumed  to  be: 


PD  =  e 


-aR2- 


where  oc  is  an  input  to  the  model  and  R  is  the  range  from  each  sonobuoy  to 
the  target.  Referring  to  Illustration  4- 13b  in  Figure  4-13  the  probability 
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Figure  4-13  PROBABILITY  DENSITY  OF  SUBMARINE'S  POSITION  WITH  RAYLEIGH 
DISTRIBUTION  IN  VELOCITY  AND  UNIFORM  HEADING:  EFFECT  OF 
PM  GRID  AS  A  FUNCTION  OF  TARGET'S  VELOCITY 


of  detection  is  assumed  to  be  of  the  fora: 


PD  -  e 


c<Pz%/v 


where  Vc  <  v 


-  C£rZ 

e 


z 


where  Vc  >  v 


where  tr  is  a  velocity  random  variable  and  Vc  is  an  input  to  the  model.  In 
effect,  this  indicates  that  the  target’s  radiated  noise  is  constant  for 
velocities  from  0  to  \/  and  that  the  radiated  noise  subsequently  increases 
linearly  as  v  >  Vc  .  In  this  example,  y  is  equal  to  9  knots. 


is  equal  to: 


For  the  case  shown  in  Illustration  4- 13a  of  Figure  4-13,  PD 


PD 


_  kPz  Vc/i5 
e 


where  tr  =  15  knots  is  the  expected  value  of  the  target's  velocity  which  is 
described  as  a  Rayleigh  distribution.  The  heading  of  the  target  is  assumed 
to  be  uniformly  distributed  between  0  and  2  ft  .  The  SPA  is  described  by  a 
joint  circular  normal  with  a  standard  deviation  of  25  n.  mi.  A  field  of  15 
sonobuoys  is  deployed  and  24  "looks"  are  taken. 

A  comparison  of  these  results  show  the  increased  capability  of 
the  sonobuoys  to  detect  the  target  submarine  as  its*  speed  increases  beyond  a 
value  tr  -  9  knots.  This  is  indicated  by  the  changes  in  probability  density 
in  the  second  quadrant  for  the  variable  PM  grid  (Illustration  4-13b)  compared 
to  that  for  the  constant  PM  grid  (Illustration  4- 13a). 


1 1 .  hffect  of  Variance  of  Variations  in  Target's  Velocity 


Figure  4-14  indicates  the  effect  of  variance  variations  in  the 
velocity  of  the  target.  The  target's  velocity  is  described  by  the  following 
probability  distribution. 
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Figure  4-14  PROBABILITY  DEMSITY  OF  SUBMARINE'S  POSITION  WITH  X  AND  Y 

VELOCITY  COMPONENTS  NORMALLY  DISTRIBUTED:  FOR  VARIANCE 
VARIATIONS  IN  VELOCITY 


In  Illustration  14a  the  variance  of  the  cr. components  of  velocity  is  64.  It 
is  equal  to  1  for  the  case  shown  in  Illustration  14b.  The  SPA  is  described  by 
a  joint  circular  normal  distribution  with  a  standard  deviation  of  25  n.  mi.  for 
both  examples  presented  in  Figure  4-14.  Also,  IS  sonobuovs  are  deployed  with 
a  total  of  24  "looks". 

2  x 

Figure  4-14  indicates  that  as  decreases  (wita  <rv* 

remaining  constant)  the  uncertainty  of  the  target's  heading  decreases.  Hence, 

when  the  velocity  of  the  target  is  described  by  a  joint  normal  and  even  when 

the  means  of  the  velocity  components  are  zero,  information  is  available  on 

the  target’s  heading  when  the  variances  of  the  velocity  components  are  unequal. 


12.  Effect  of  Biasing  Target  Heading 


Figure  4-15  presents  results  obtained  for  the  case  in  which 
target  velocity  is  described  by  the  following  probability  density: 


f  (ir£,Vy) 


_L _ „  f 

ZTr  P  [  2  \  1*  +  1¥V  / } 


The  initial  SPA  is  described  by  a  joint  circular  normal  with  a  standard 
deviation  of  25  n.  mi.  A  field  of  15  sonobuoys  are  deployed  within  the  SPA. 
A  total  of  24  "looks"  are  taken  during  the  elapsed  time  interval  between 
3  to  9  hours.  Each  sonobuoy  is  monitored  on  every  other  "look". 


The  effect  of  inserting  a  mean  value  of  12  knots  to  the 
component  of  target  velocity  is  to  bias  the  heading  of  the  target.  The  joint 
probability  density  of  the  target's  velocity  and  heading  is  given  by: 

.  ir  -(vz+Uz)/2er1  cos  &/cr i 

3<K0>‘I^e 

(This  equation  is  presented  in  Reference  3,*  page  196)  where  cr1  =  144  and 
u  A  =  12.  Integrating  out  &  results  in  the  following  velocity  distribution: 

**/>{“  (rr^)/*<r*}  j  (Mjl) 

‘Reference  3:  Probability  Random  Variables  on  Stochastic  Processes  by 
Athanasios  Papoulis  Published  1965  by  McGraw-Hill. 
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Figure  4-15  PROBABILITY  DENSITY  OF  SUBMARINE'S  POSITION  WITH  X  AND  Y 

VELOCITY  COMPONENTS  NORMALLY  DISTRIBUTED:  EFFECT  OF  BIASED 


Curves  of  ^(za)  and  the  heading  distribution  y(£)  are  presented 

u.  7 

in  Reference  2  on  page  562  and  365  for  a  number  of  values  of  /  =  —2—  . 

2crz 

The  results  shown  in  Figure  4-15  indicate  that  the  effect  of 
biasing  submarine  heading  is  to  shift  the  high  density  area  of  probable 
target  positions  to  a  sector  located  outside  of  the  ring  of  deployed  sonobuoys. 

13.  Accumulated  Probability  Density  of  Target  Position  with  Search 
Time 

A  limited  effort  has  been  conducted  toward  development  of  a 
program  which  computes  the  probability  that  the  target  will  pass  through  :» 

Ax  by  Ay  grid  during  the  search  time  from  Tf  to  Tz  .  Some  sample 
results  of  this  effort  are  presented  in  Figures  4-17  and  4-18.  For  Figures  4-17 
and  4-18  the  velocity  of  the  target  is  assumed  to  be  13  knots.  The  heading 
of  the  target  is  described  by  the  following  distribution: 

p(e>T~ “*?{“**} 

The  accumulated  time  is  the  interval  between  3  to  9  hours.  Figure  4-17 
presents  results  for  the  case  in  which  no  sonobuoys  are  deployed.  Figure  4-18 
presents  results  for  the  case  in  which  one  sonobuoy  is  deployed. 

The  example  presented  in  Figure  4-18  is  the  same  as  that 
shown  in  Figure  4-17  with  the  exception  that  target  heading  is  assumed  to  be 
uniformly  distributed  over  the  interval  from  0  -  -g  .  This  angle  is 
measured  in  the  clockwise  direction. 
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SYMBOLS  INDICATE  VARIT'JS  PROBABILITY  DENSITY  LEVELS 
DARKER  AREAS  INDICATE  HIGHER  PROBABILITY  DENSITY  AREAS 
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O  datum 


-J 

SUBMARINE  HEADING:  p{-&)=  -  exp(cos€>-\ 

*/7l0(1)  f 

SUBMARINE  VELOCITY:  13  KNOTS 

TIME  LATE:  3  HOURS 

ACCUMULATED  TIME:  3-9  HOURS 

GRID  SIZE:  2.5  x  2.5  N.MI. 

NUMBER  OF 

SONOBUOYS:  1  ,  NUMBER  OF  "LOOKS":  24 

Figure  4-17  ACCUMULATED  PROBABILITY  DENSITY  OF  SUBMARINE'S  POSITION  WITH 
KNOWN  CONSTANT  VELOCITY  AND  RANDOM  HEADING:  ONE  SONOBUOY 
DEPLOYED 
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SYMBOLS  INDICATE  VARIOUS  PROBABILITY  DENSITY  LEVELS 
DARKER  AREAS  INDICATE  HIGHER  PROBABILITY  DENSITY  AREAS 
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O  DATUM 


SUBMARINE  HEADING:  UNIFORM  DISTRIBUTION  0-  272 

SUBMARINE  VELOCITY:  13  KNOTS 

TIME  LATE:  3  HOURS 

ACCUMULATED  TIME:  3-9  HOURS 

GRID  SIZE:  2.5  x  2.5  N.MI. 

Figure  4-18  ACCUMULATED  PROBABILITY  DENSITY  OF  SUBMARINE'S  POSITION  WITH 
KNOWN  CONSTANT  VELOCITY  AND  UNIFORM  HEADING  BETWEEN  0  AND 
7T  /2:  NO  SONOBUOYS  DEPLOYED 
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V. 


APPENDIX  A:  CAL  FLYING  SPOT  SCANNER 


The  CAL  Fiying  Spot  Scanner  was  developed  as  a  flexible  research 
facility  for  investigation  of  pattern  recognition/ image  analysis  problems. 

It  has  been  interfaced  with  various  computers  at  the  Computer  Center  as  these 
have  been  updated  from  generation  to  generation  and  is  currently  capable  of 
on  line  operation  with  either  an  IBM  570-165  or  a  PDP-9  computer.  General 
software  utility  programs  have  been  developed  for  reading  and  writing  (computer 
input/output)  operations.  Recently,  a  new  utility  program  has  been  developed 
to  adapt  the  standard  IBM  570  plotting  routines,  provided  for  use  with  electro- 
optical  and  electro  mechanical  plotters  at  the  Center,  to  permit  direct  control 
of  the  Flying  Spot  Scanner  as  a  plotter. 

Plotting  instructions  to  the  Flying  Spot  Scanner  can  now  be  oirectiv 
executed  by  FORTRAN  and  PL/I  programs.  Included  is  the  capability  to  generate 
the  Standard  character  and  symbol  set.  This  feature  together  with  the  variety 
of  photographic  recording  devices  incorporated  in  the  scanner  has  greatly  expanded 
the  utility  of  the  system. 

The  Flying  Spot  Scanner  features  a  cathode  ray  tube  controlled  by  a 
precision  magnetic  deflection  system.  Nominal  spot  size  is  .0018"  with  1024 
discrete  addresses  per  line  horizontally  and  1024  lines  vertically.  No 
programming  restriction  is  placed  on  successive  spot  address  location,  a  feature 
that  facilitates  operation  in  the  plotting  mode  as  well  as  in  more  complex 
interactive  pattern  reading  modes.  Two  complete  cathode  ray  tubes  and  deflection 
systems  are  provided  which  can  be  connected  alternatively  to  common  logic  control 
and  power  supplies.  One  is  configured  to  read  opaque  copy  and  is  used  only  in  the 
read  mode.  The  second  CRT  is  designed  to  read  transparencies  and  to  write  on 
photographic  material.  A  plate  holder  is  provided  to  mount  transparencies  with 
a  useful  scan  area  of  5"  x  3".  A  film  transport  which  features  pin  regis¬ 
tration  to  an  accuracy  of  .001"  under  computer  control  is  provided  for  reading 
70mm  and  35mm  film.  For  recording  purposes  a  variety  of  cameras  and  cut  film 
holders  as  well  as  Polaroid  backs  are  provided.  A  35mm  and  a  16mm  movie  camera 
with  single  frame  advance  control  by  computer  may  be  used  alternatively.  Pro¬ 
vision  is  made  for  introducing  color  filters  in  the  optical  path  thus  allowing 
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color  separation  in  the  read  node  and  color  writing  with  appropriate  color 
negative  naterial.  A  large  screen  nonitor  CRT  display  is  provided  to  observe 
the  action  of  the  scanner  in  the  reading  node  or  to  display  in  parallel  the 
information  being  written  on  the  photographic  negative. 


A  VIEW  OF  THE  FLYING  SPOT  SCANNER  SHOWING  THE 
CRT  AND  OPTICAL  BENCH  IN  OPENING.  TOP  LEFT.  AND 
THE  LARGF  SCREEN  MONITOR  C-M  RICH:. 


ILLUSTRATIONS  OF  THREE  DIMENSIONAL  DATA  PLOTS  BY 
370/FLYIN6  SPOT  SCANNER  SYSTEM 


LINE  PLOT  (35  mm) 


